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Abstract: A pedagogical introduction to color superconductiv- 
ity in the weak coupling limit is given. The focus is on the ba- 
sic tools of thermal field theory necessary to compute observ- 
ables of color superconductivity. The rich symmetry structure 
and symmetry breaking patterns are analyzed on the basis of 
the Anderson-Higgs mechanism. Some techniques can also be 
applied for computing neutrino processes in compact stars. 
As an example, we show how to obtain the neutrino emis- 
sivity for Urea processes in neutron stars by computing the 
polarization tensor of the W-boson. We also illustrate how a 
spin-1 color superconducting phase generates an anisotropic 
neutrino emissions in compact stars. 
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1 What is color superconductivity? 

Nowadays, superconductivity in a wide range of materials and superfluidity 
of Helium-3 can be generated with practical equipments in laboratories. It 
is well known that the underlying mechanism of all these 'super' phenom- 
ena is the so-called Cooper paring of fermions (electrons in superconductors, 
fermionic 3 He atoms in Helium superfluids) near the Fermi surface of the cor- 
responding fermions. These phenomena are all of electromagnetic nature and 
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at atomic level. In the regime of strong interaction where quarks and gluons 
are elementary particles described by quantum chromodynamics (QCD), a 
similar phenomena, the so called color superconductivity was first proposed in 
cold dense quark matter by Barrois [l], Bailin and Love p], and further devel- 
oped by many others (for reviews, see |3l HI El El d El El llffl ED C21 C31 E]). 
However color superconducting quark matter cannot be produced in any 
laboratory in the foreseeable future. The reason is that this exotic phase 
of matter requires extremely high baryonic densities and relatively low tem- 
peratures. In nature, such conditions are only realized in the cores of cold 
compact stars, i.e. in relatively old remnants of supernova explosions. 

In a color superconductor, quarks form Cooper pairs similarly to electrons 
in a normal superconductor or to fermionic 3 He atoms in a Helium superfluid. 
Quarks carry color charges and interact via gluon exchange (the gluon is the 
QCD gauge boson). At very high baryonic densities the quark chemical po- 
tential fj, is much larger than the QCD scale Aqcd ~ 200 MeV, the running 
coupling constant as(fi) is thus small, because QCD is an asymptotically 
free theory for quarks and gluons. This was first shown by Gross, Politzer 
and Wilczek in 1973 pi E] and awarded the Nobel Prize in 2004 [T7HT3] . 
In this case, the interaction between quarks is dominated by one-gluon ex- 
change. As will be demonstrated explicitly at the end of this section there is 
an attractive channel in one-gluon exchange providing the binding agent for 
the quark Cooper pair. This is in contrast to normal superconductors where 
an electron Cooper pair is formed indirectly via electron-lattice interaction 
(direct electron-electron interaction is repulsive). In this sense color super- 
conductivity is simpler than the normal one. In hadronic matter quarks are 
confined in color-neutral hadrons. Only at sufficiently high baryonic densi- 
ties when these hadrons overlap and eventually free the quarks out of their 
hadronic domains, the quarks can move in much larger space-time volumes 
than the size of hadrons. This phenomenon is called deconfinement. It is ex- 
pected to occur in the core of some compact stellar objects as neutron stars, 
where the baryonic densities are a few times larger than the nuclear matter 
saturation density po ~ 0.14 fm -3 , which corresponds to a quark chemical 
potential y. ~ 400 - 500 MeV. 

Since quarks are fermions, the Pauli exclusion principle does not allow 
two quarks to occupy the same quantum state. At zero temperature, they 
fill up states up to the so-called Fermi momentum kp (or the Fermi surface) 
above which all states are empty. Due to Pauli blocking, only quarks near the 
Fermi surface are able to interact and exchange momenta. The exchanged 
momenta in quark-quark scattering are of order kp. In the ultra-relativistic 
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limit, we have Ep = J kp + m 2 « kp ~ fj,. We know that the quark density 

p ~ kp ~ /i^, the high baryonic density means large Fermi momentum kp or 
chemical potential //. We will consider asymptotically large quark chemical 
potentials, /x ^> Aqcd, the so-called weak coupling approach, which means 
high density and large mementum transfer. The weak coupling approach 
enables one to study the phenomena in a rigorous or well-controlled way. 
For reviews of the weak coupling appoach, see, for example, [5[ [13]. The 
physical predictions of such an approach certainly have to be interpreted 
very carefully, since for realistic \i one has a s (fi) ~ 1. But aided by the 
resummation methods, the perturbative analysis can be extrapolated even 
to realistic densities such as those in the cores of neutron stars [13J. Support 
also comes from the NJL model [El ED], a simple effective theory of QCD 
without gluons but with a local four-quark interaction. It predicts color 
superconductivity also at moderate baryonic densities with gaps of order 
of (j) ~ 100 MeV, which matches the gap value extropolated from the weak 
coupling approach. Recently there are many developments in the NJL model 
[2T1 1221 1231 [24] and in the random matrix model [25, 26J, a more rigorous 
effective model, to describe the color superconducting phase diagram. 

So at high baryonic densities where the strong coupling constant is small, 
the single gluon exchange dominates the quark-quark scattering, as shown 
in Fig. ([T]), whose amplitude is proportional to 

N '~ l iV c + l 

E T a' T "' = Jy— - Vi'j) 

+ -^- (Sii'fyj' + Sij'Si'j) ■ (1) 

Here, N c = 3 is the number of colors and T a are generators of the SU (3) c 
gauge group. Here T a = X a /2 where A a are the Gell-Mann matrices. The 
indices i, j are the fundamental colors of two quarks in the incoming chan- 
nel, while i',f their resprective colors in the outgoing channel, see Fig. ([T|. 
Interchanging two color indices in the incoming or outgoing channel changes 
the sign of the first term while the second term remains intact. The minus 
sign in front of the asymmetric term indicates that this channel is attractive, 
as it is similar to the Coulomb interaction between a negative electric charge 
and a positive one (the product of two charges is negative). This finding is 
crucial, as due to Cooper's theorem any arbitrarily weak attractive interac- 
tion will destabilize the Fermi surface in favor of the formation of Cooper 
pairs. Since the pairs are of bosonic nature they will condense at sufficiently 
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low temperatures into a Bose-Einstein condensate (in a more rigorous way, 
Cooper pairs are Bose-Einstein condensate in strong coupling limit). 

In group theory, Eq. (Q]) corresponds to deducing the direct product of 
two triplets into one color antitriplet (anti-symmetric) and one color sextet 
(symmetric) 

[3] c ®[3] c = m®m s c . (2) 

Now we illustrate how two triplets makes an anta-triplet. Denote ipi the 
quark field following the SU(3) transformation: 

4>'i = Qijipj, (3) 
An anti-triplet can be written by ipk = eijktpiipj, 

- dk'ktijk'gii'gjj'ipi'^j' 
= gj k det(g)e i > j > l ip i >'ip j > 

-- vi, (4) 

where we have used tijk'gu'gjj'gk'i = det(g)ei'j'i and det(g) = 1. We see that 
the anti-symmetric diquark field ipk transforms like an anti-quark. 

From the anti-symmetry of the attractive channel it follows that the two 
quarks in a Cooper pair must carry different colors. In addition to colors, 
quarks also come along with spin and flavor quantum numbers. Since the 
total wave function of a Cooper pair has to be anti-symmetric under the 
exchange of two quarks, the combined spin-flavor part has to be symmetric. 
This requirement will be used in the following section to classify various 
possible color superconducting phases. 

2 Structure of diquark condensate 

In this section we discuss the structure of diquark condensate [HI [27]. A 
diquark condensate is defined as an expectation value 

M = (? C 0V> , (5) 



where g is a SU (3) element, 
which transforms as 

^'k ^ 
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Figure 1: One gluon exchange 




where ip is the quark field with all Dirac, color and flavor indices, ip c is 
its charge conjugate field defined by ip c = ip T C with the charge conjugate 
operator C = z7 2 7o- The operator O has Dirac, color and flavor parts: 

O = O Dirac ® O color ® O fl avor . (6) 

The charge conjugate operator can be determined by transforming the Dirac 
equation of a fermion to that of an anti-fermion, 

{Piil^ - m)V> = 
^ T {P^ T -m) =0 
-»• i?C{ Pll C- l ^ T C - m) = 

- f_ c (p»C- 1 1 » T C-m) = 

- Mp^ + m) = 0. (7) 

One can verify that C = i7 2 7o satisfies C~ l ^ T C = -7^ or C^ T C = 7^ 
using 7^ T = 7 , — 7 1 , 7 2 , —7 s for ji = 0, 1, 2, 3. Exploiting the anti-symmetric 
property of the fermionic quark fields, one finds 

- -{^{CO)l^) = -{^ T {CO) T ^). (8) 
So the operator CO obeys 

(CO) T = -CO. (9) 

This means that the operator CO is anti-symmetric under transposition. 
We can choose all three parts of CO, i.e. Dirac, color and flavor parts, anti- 
symmetric or two of them symmetric and the third one anti-symmetric. In 
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Table 1: Symmetric properties of O under transposition. In the row of Dirac 
part, 'S', 'P', 'V, 'A', 'T' means that the operator CO is a scalar, pseudo- 
scalar, vector, axial-vector and tensor, respectively. In the U(2) row, <7j are 
Pauli matrices, while in the U(3) row, Aj are Gell-Mann matrices. 





Anti-symmetric 


Symmetric 


Dirac 
(combined with C) 


75 1 7 M 75 
S P V 


A T 


U(2) 


singlet [1] 


1, CT3 
triplet [3] 


U(3) 


A2, A5, A7 
anti-triplet [3] 


1, Ai, A3, A4, A6, As 
sextet [6] 



Tab. {T]) a variety of choices to build up the operator O is given together 
with the corresponding symmetric properties under transposition. 

Now we show how to determine the parity of the condensate. Suppose 
under parity transformation V, the field ip{t, x) changes as follows 

P^(t,x)V = 7770^,-x) 

7ty(t,x)P = r ? >(t,-x) 7o , (10) 
where 77 is a complex number. Then the conjugate fields transform as 
Vtp c (t,yi)V = rjcjoipcit, -x) 

P$ c {t,yL)V = r£fo(t,-x) 7o . (11) 
So the condensate transforms as 

P^ c (t,x)0^(i,x)P = v*M c {t, -x) 7 o0 7 oV(t, -x) 

= -^ c (t 5 -x) 7o 7 o^(t,-x), (12) 

where use was made of the fact that the parity of the conjugate particle is 
opposite to that of the particle, i.e. r]* = —77. Then O = 75, 1, 7 M 7 5, 7 At , a**" 
correspond to a scalar, pseudo-scalar, vector, pseudo-vector and tensor, re- 
spectively, as shown in the second row of Tab. (pQ. For example, O = 7^75 
gives a vector, since 

P^ c (t,x)7%s^(i,x)7> = ^ c (t,-x)7°75^,-x), 
^ c (t,x)y 75 V'(i,x)P = -^ c (t,-x) 7 i 75^(i,-x) ; (13) 
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Table 2: Structure of diquark condensates M, the Dirac part is refered to 
that of the operator O (not that of CO). J p is the spin-parity. Note that 
the color part is always in anti-triplet for all condensates. 



pairing 


Dirac part (J p ) 


flavor multiplet (number of flavor) 


spin-0 


75(0 + ),l (0-) 


1 (N f = 2),3(N f = 3) 


spin-1 


r (I") 


1 (N f = 1), 3 (N f = 2), 6 (N f = 3) 


spin-1 


7^75 (1+) 


1(JV> = 2),3(JV> = 3) 



which transforms as a vector. 

Prom spin-0 pairing (see, e.g. [IH [20j [27l EH]), we require that the color 
part is anti-symmetric with respect to exchanging two quarks, i.e. the color 
part is in the anti-triplet channel, and that the Dirac part must be a Lorentz 
scalar or pseudo-scalar, i.e. the Dirac part of O must be 75 or 1 which is 
anti-symmetric combined with C . Therefore the flavor part should be anti- 
symmetric too, i.e. it must be a flavor singlet in the two-flavor case or a flavor 
anti-triplet in the three- flavor case. The condensate is then a J p = + (spin- 
parity) or J p = 0~ bound state of two quarks depending on whether the 
Dirac part is 75 or 1. For the spin-1 pairing (see, e.g. [29j [30j EU EH 
EH EH EH El]), the color part must be anti-triplet (anti-symmetric) and the 
Dirac part can be a vector (symmetric) or an axial- vector (anti-symmetric), 
so the flavor part can be symmetric or anti-symmetric respectively. There 
is another example for the spin-1 pairing in the single-color case where the 
color part is symmetric and flavor one anti-symmetric |38| . If the single 
flavor pairing occurs (the flavor part is symmetric), the Dirac part must be 
7 M . Therefore a spin-1 Cooper pair is a J p = l - state. Tab. J2]) summarizes 
the structure of diquark condensates in spin-0 and spin-1 pairings. 

Generally the gap of a spin-0 pairing is much larger than that of a spin-1 
pairing. In the ideal case there is more energy benefit for the spin-0 pairing 
compared to the spin-1 pairing, so the spin-0 pairing is favorable. But in 
the real world, there are many factors to disfavor the spin-0 pairing. For 
example, when electric neutrality is taken into account, there is a difference 
between the chemical potential of the strange quark and that of light quarks 
u and d due to the large strange quark mass. Additionally, the chemical 
potentials of u and d quarks differ due to /3-equilibrium. If the differences of 
the chemical potentials are large enough, the spin-0 pairing, e.g. CFL phase, 
is not favorable [39j [40j [4TJ [42] , and the spin-1 pairing or LOFF pairing [43] 
might be the true ground state [44], [45] . Another possibility to kill the spin-0 
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Figure 2: The ferromagnetism as an example of spontaneous symmetry 
breaking. 




pairing is the presence of a strong magnetic field which favors the symmetric 
spin wave function [32] . 

3 Spontaneous symmetry breaking 

A symmetry of the Lagrangian is said to be spontaneously broken if the 
ground state or the vacuum of the system is not invariant under the oper- 
ations of that symmetry. We call this phenomenon spontaneous symmetry 
breaking. 

A well-known example is the ferromagnetism arising from the spin-spin 
coupling, see Fig. J2|). The Lagrangian describing the ferromagnetism is 
invariant with respect to a SO(3) rotation. Above the transition temper- 
ature T c , the vacuum state of the spin system is invariant with respect to 
50(3) rotation since the spin orientations are totally random and any rota- 
tion does not change the spin state macroscopically. Below T c the vacuum 
spontaneously chooses one magnetization direction, so the groud state is no 
longer invariant under 50(3) rotation. The 50(3) symmetry of the vacuum 
is spontaneously broken to 50(2) which describes the rotational symmetry 
around the total spin direction. We thus write this symmetry breaking pat- 
tern as G = SO (3) — ► H = 50(2), where G is the total symmetry and H 
the residual symmetry. 

In the second example of spontaneous symmetry breaking the concept of 
the locking state is explained. Suppose we have a system composed of red and 
blue spins of U{1) symmetry, see Fig. j3|). Above a transition temperature 
T c , the red and blue spin orientations are totally random. Below T c , there 
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Figure 3: The locking state as an example of spontaneous symmetry break- 
ing. 




are four possibilities of breaking the total symmetry spontaneously: 

G = U{l) red U{l) blue 

G = U(l) red ® U{l) b lue 

G = U(l) red ®U(l) blue 
G = U(l) red ® U(l) Uue 



H — U(l)blue 

H = U(l) red 
H = 1 

H = U(l) red+ blue 



(14) 



In the first one the red spin symmetry is broken while the blue one remains. 
In the second one the blue spin symmetry is broken but the red one remains. 
In the third breaking pattern, both the red and blue spin symmetries are 
spontaneously broken. The fourth pattern is special in that both the red 
and blue spins point in any directions but their relative angle is fixed. This 
symmetry breaking pattern is called locking in the sense that the state is 
invariant under the joint rotation of the red and blue spins with the same 
rotational angle. The locking, for example, is the underlying mechanism for 
various phases in the superfluid Helium-3 [46J. 

The well-known Nambu-Goldstone theorem [EJUS] states that the spon- 
taneous breaking of a continuous global symmetry implies a massless spin 
zero or scalar particle, and that each broken generator of the symmetry 
group gives rise to one such particle. Consider the complex scalar field as an 
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example. The potential can be written as 

v{\fa 2 ) = a\fa 2 + b\fa\ (15) 

where a ~ a'(T — T c ) and a',b > 0, see Fig. The above potential is 

invariant under £7(1) transformation, fa = e ld (f). Above T c the potential ex- 
hibits a minimum at <p = which is intact under £7(1) transformation. Below 
T c the potential is minimal along the circle with radius \fa = fa, because 
the coefficient a has different signs below and above T c . The symmetry of 
£7(1) is spontaneously broken because the vacuum is changed after any £7(1) 
transformation fa = e l8 (j) = e td fa. One sees that below T c there is a massless 
mode (the so-called Goldstone mode) along the circle of minima. 

A general proof of the Goldstone theorem can be given as follows. Sup- 
pose a given potential is invariant under the transformation of a continuous 
group: 

5<S> = i9 a C a $, (16) 

where $ = [<f>\, ■ ■ ■ , <ft nG ] T with fa being the real scalar field component and 
«G the number of real scalar fields. C a are operators of the corresponding 
Lie algebra of the group. The infinitesimal change of the potential due to 
5$ should vanish 

6V (®) = YT 5 ^ = i ^r ea ^ ca )^ = ^ (17) 



Since 9 a are arbitrary we have 

SV 



(C a )ii^ = (U 



for a = 1, • • • ,na- Taking an additional derivative with respect to fa, we 
obtain 

5 2 V 5V 

-(C a hfa + —{C a ) ik = 0. (19) 



'>kO<Pi 



Evaluating the above equation at the vacuum $ = $o, the second term 
vanishes, and we have 

MUC a h<t>0j = 0, (20) 



where the mass matrix is defined by = A J, 



. If the ground state 

<I>='I>0 

is invariant under a sub-group H of G, for each generator belonging to this 
sub-group H, we have 



(C^ijfaj = (21) 
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Figure 4: Potential for complex scalar field 



v(\4>\ 2 ) 



A 




A 



Re(0) 



T > T c 



T <T C 



c 



for a = 1, • • • , n# . For the remaining no — nn generators, we have 



So Eq. (f20l) shows that there are uq — nn zero eigenvalues for the mass 
matrix which means no — nn massless Goldstone bosons. 

Note that the above argument holds for a Lorentz invariant system. In 
the context of dense quark matter where the Lorentz invariance is lost, the 
Nambu-Goldstone bosons have many fine structures, see, e.g. Ref. [IH [50l 



4 Scalar QED as an example of superconductor 

In the above section we discussed the spontaneous breaking of a global sym- 
metry. In this section we will investigate the spontaneous breaking of a local 
gauge symmetry. We will show that it is also a good example for supercon- 
ductivity. To this end, we consider a theory of a complex scalar and a U(l) 
gauge field where they couple to each other in a minimal way, the so-called 
scalar quantum electromagnetism or scalar QED. The Lagrangian is 



(22) 



ED. 



c = {d^d^-vm 2 )-\f^ 



(23) 
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where the covariant derivative is defined by = <9 M + ieA^ and the guage 
field stength tensor by F^ v = 8^A U — 8 u A fM . The potential is given by 
^(l^l 2 ) = -^(l^l 2 ~~ 0o)- The Lagrangian £ is invariant under a U(l) local 
gauge transformation, 

<f>' = e ied(x U, 

A'^ = A^-8^9. (24) 

One can verify that under the above transformation D'^cp' = e ie6< ^ D^cj). 
Therefore the term (D^y D^cf) is invariant. 

The Lagrangian equation for the gauge field reads 

8^ F^ u = j u , 

j v = ie(4>*d v (t>-(t>d v (t>*) -2e 2 \4>\ 2 A u . (25) 

The current satisfies the conservation law <9 M j M = due to 8^8 V F^ V = 0. 
The Lagrangian equation for the scalar field reads 

{d^ + ieA^ + ieA^ = 2A0(0§ - |0| 2 ). (26) 

In the covariant gauge d^A^ = and choosing the real scalar field (unitary 
gauge) (p* = (j) = p, the equation for the gauge field becomes 

(d^ + 2e 2 p 2 )A u = 0. (27) 

Suppose p is very close to the vacuum value (fio, we can expand p = (fro + r/, 
the linearized equation of motion for the scalar field become 

(d^ + 4\cf> 2 )r} = 0, (28) 

where we have dropped higher order terms like A^d^r] and A^A^ etc.. 



5 Invariant subgroup for order parameter 

Consider a symmetry group in a theory 

G = Gi®G 2l (29) 

where G\ and G2 are Lie groups (a generalization to more than two groups 
is straightforward). The group elements can be expressed in terms of the 
generators of the respective Lie algebras 

9 = 9i®92 
gi = e iC ^eG 1 

g 2 = e iC ^EG 2 , (30) 
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where C" (a = 1 • • • n%) and Cf (a = 1 • • • 77-2) are the generators of the Lie 
algebras G\ and G 2 , respectively. Of and Q% are rotational angles associated 
with the operators Cf and Cf respectively. 

Now we investigate how the order parameter changes under the transfor- 
mation of the group G. After choosing one representation for the group G, 
the order parameter M. can be written on its basis 

M rsj Ayeii <g> e 2 j, (31) 

where the coefficient Ajj is the gap matrix, and ey and e 2 j are basis vectors 
for Gi and G 2 representations, respectively. After transformation, the order 
parameter becomes 

M' = gx® giMgl <g> gf 

~ AijgxeugJ <g> g 2 e 2 jgl '• (32) 

For infinitesimal transformation 

51 « 1 + *C?0? 

52 « l + iC 2 % a , (33) 

Eq. JHH) becomes 

M' = A^ + ^A, J (C7 1 a e l4 + e lj C 1 aT )»e 2j 

+t^Ai i ei i (C7 2 a e 2j + e 2j C 2 aT ) (34) 

To find the invariant subgroup for the order parameter, we require M' = M 
and then 

ie a 1 A lJ (C'?e li + e li Cf)®e 2j 

+i^Ayeii (8) (C 2 e 2i + e 2i C7 2 aT ) = (35) 

From this equation we can find the subset of operators Cf 2 which make up 
the sub-algebra of the residual symmetry. 

For the SU(3) group, we have a simpler version of Eq. (135!) . We know 
that the color part of the order parameter is the anti-triplet. So the basis of 
the color part now should be Jj, the color anti-symmetric matrices defined by 
(Ji)jk = iejik = —itijk- Now we look at g\J%g\ , where g\ is a SU(3) group 
element in fundamental representation (a 3 by 3 unitary matrix), which 
satisfies gig\ = 1 and det(gi) = 1. Using these properties we can simplify 
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giJigJ as 

9lji[giJi9l]mn = 9ljigiml{Ji)l P ginp 



— V€ilp9lji9lml9lnp 

= -ie jmn det(gi). (36) 
Multiplying (g\)kj on both sides, we get 

[9l Jk9l ]mn = ~itjmn{9l)kj 

= (9l)kj(Jj)mn, (37) 

or in operator form 

9iJk9i = (g\)kjJj- (38) 

In the case of the SU(2) group, for example, the number of flavors is 2, 
we will encounter gi<J2gJ , where g\ is a SU(2) group element in fundamental 
representation (a 2 by 2 unitary matrix). Using g\g\ = 1 and det(gi) = 1, 
we have 

[gi^gTjmn = g\ml{02)l P g\n V 

= -iei p gimiginp 
= -ie m „det(#i) 

= (^)mn (39) 

Therefore we obtain gi<J2gJ = o-i- 

We can also simplify g2^jg2 f° r SO(3)j. Here Kj with j = 1, 2, 3 are the 
basis for the spin triplet. Because 52 is a group element of SO (3), we have 

rp 1 

92 = 92 ■ Then we obtain 

92^j92 = 92KJ92 1 

~ Kj ~\~ %0<2 K^nKj ^^2 ^j^n 
= Kj + 102 i^mjl^l 
= Kj — 1&2 {Jm)jl K l 
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5.1 Symmetry patterns: 2SC phase 

Our first example for the spin-0 pairing is the 2SC phase with u and d quarks. 
For simplicity we consider the positive parity channel, for the complete case, 
see, e.g., Eq. (9) in Ref. [14]. The order parameter for this color- flavor 
coupling phase can be written in the following form 

M = JiAijIjjs, (41) 

where Jj is the color anti-symmetric matrices given by (J{)jk = itjik = 
—ieijk- We can also express Jj in terms of Gell-Mann matrices J\ = A7, 
J 2 = — A5 and J3 = A2. I{ are the flavor anti-symmetric matrices defined 
in the same way as J«. The coefficient of the order parameter Ay for 2SC 
is Ay = A5js5i3. Incorporating this expression for A^- into Eq. (JHJ, we 
obtain 

M 2SC = AJ 3 fj 2 7 5 , (42) 

where (a 2 )ij = {h)ij = -i^ij with i,j = 1,2. 
The symmetry group for the 2SC phase is 

G = Gi®G 2 ®G 3 , (43) 

where G\ = SU(3) C is the color group, G 2 = SU(2)f the flavor group, and 
G3 = U(1)b the baryon number group. Note that the electromagnetism 
group U(l) em is a subgroup of SU(2)f and U(1)b, since the electric charge 
satisfies 

Q = diag(Q u , Q d ) = diag(2/3, -1/3) = T 3 + B/2, (44) 

where B = 1/3 is the baryon number. 

First we want to find the invariant subgroup of SU (3) c and SU (2)j, so we 
can safely drop G3 for simplification of the notation. The symmetry group 
now is 

G = Gi®G 2 = SU(3) c ®SU(2) f . (45) 
The elements of SU(3) C and SU(2)f are 

9 = 9i ® 92 

92 = e^eSU(2) f , (46) 
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where a a with a = 1,2,3 are Pauli matrices. According to Eq. (|32l) . we 
obtain 



M' = gi<S> g 2 MgJ <g> g 2 

= A7551 J 3 gJ ® g 2 a 2 gl 

= A^ 5 (g\) 3i Ji ®<t 2 (47) 

where we have used Eq. ([38]) and (|39l) . By expanding g± and 52 to the linear 
term of Of and Q\ respectively, and requiring M.' = M., we have 

et{T a ) M J i ®G 2 = 0, (48) 

For a = 1,2,3 Eq. (14"51) holds automatically, because elements on the third 
row of T\,T 2 ,T% are zero. This means that there is freedom to choose any 
values of 9f with a = 1, 2, 3. For a = 4,5, 6, 7, 8, we have 

WTjuJi = {e\ + i6\)T 7 -{el + iel)T b -^=elT 2 (49) 

where we see that 9^{T a ) 3 iJi is not zero. Therefore to make Eq. (|48l) hold, 
0" must vanish. We see that the orginal SU(3) C symmetry is spontaneously 
broken to SU(2) C spanned by operators T\, T 2 and T3 as follows 

SU(3) c ®SU(2)f -» SU(2) c ®SU(2) f . (50) 

In order to get the invariant U(l) group for the order parameter, we 
consider the transfromation of U(l) em for electromagnetism in the flavor 
sector. The elements of SU(3) C and U(l) em are 

g = g\®g 2 

91 = e^eSU(3) c 

g 2 = £U(l) em , (51) 



-01(T a ) 3i Ji ® a 2 + J 3 ® a 2 j = 0. (52) 



Eq. (1351) becomes 



where we have used 



-ie fcZ det[e^] 

(Q,+Q d )e = a2e ie/s (53) 
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Table 3: Q charge in 2SC and CFL phases. The colors are labelled by 1,2 
and 3. 



2SC 


1 


2 


3 


u 


1/2 


1/2 


1 


d 


-1/2 


-1/2 






CFL 


1 


2 


3 


u 





1 


1 


d 


-1 








s 


-1 









The above equation becomes 

[^ e \-ie\)T 7 + {et + t9\)n 

2 

We obtain 

na 



2 

® o-2 + -0 2 T 2 ®<T 2 = 0. (54) 
o 



Ul _ 0, a = 4, 5, 6, 7 
^ = ~^f 2 . (55) 

With this relation between Of and 02 we can combine the group elements as 

j*** toe*** = (56) 
with the new U{1) charge given by 

Q = 1 <g> Q - -^T 8 (8) 1, (57) 
v3 

where Tg = ^=diag(l, 1, —2). So the original symmetry is broken as 

SU(3) C ® C/(l) em -> 5C/(2) C ® E/(l). (58) 

The Q charge in the 2SC phase for u and d quarks with different colors are 
listed in Tab. J3]). We see that the condensate made of (uid2,U2d\) is Q 
neutral because Q ccm d = Q{u\) + Q{d 2 ) = Q(u 2 ) + Q{d\) = 0. 

5.2 Symmetry pattern: CFL phase 

Our second example for the spin-0 pairing CFL phase. For simplicity we 
consider the positive parity channel, for the full case, see, e.g., Eq. (6) in 
Ref. [15 • The order parameter can be written as 

M = JiAijljis = AJ ■ I 75 , (59) 
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with Aij = A5ij. The symmetry group is 

G = GtxG 2 xG 3 = SU{3) C x SU(3) f x U(1) B - (60) 

Now we try to find the invariant subgroup for the order parameter. We will 
see U(1)b plays no role, so we only G\ and G2 in the following discussion. 
Then the group elements are 

9 = 91 ® 92 

91 = e^€SU(3) c 

92 = e'^eSP(%. (61) 
Eq. (|35l) for the CFL phase becomes 

+0%Ji (Tjijlj = 0. (62) 
We use following formula to simplify the above equation 



P(Ta)ijJj = (-9 1 + i0 2 )T 5 + (0 A - i9 5 )T 2 

v3 

a (T a )2jJj = {e l + i6 2 )T 7 + {9 & -id 7 )T 2 



1 

e^r^j, = (e 4 + ie 5 )T 7 -(e 6 + ie 7 )n 

~^0 8 T 2 (63) 

Then we obtain following conditions from Eq. (|62l) 

91 = VaOZ (64) 

where r] a = 1 for a = 2, 5, 7 and rj a = —1 for a = 1,3, 4, 6, 8. We see that in 
order to make the order parameter invariant, we have to lock the rotaional 
angles for SU(3) C and SU(3)f. The orginal symmetry SU(3) C x SU(3)f is 
broken to SU(3) c+ j. 

To find the new U(l) charge, we consider the transfromation of U(l) em 
in the flavor sector with the generator of U(l) em given by 

Q = I Q d I = I -1/3 I . (65) 
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Note that U(l) em is a subgroup of SU(3)f since Q = T3 + Tg/y^. We use 

= -ie ikl det[e iQe ] 

= j ie *(Q»+Qd+Q.)0 = j. ( 66 ) 

to obtain 

e iQ0j. e iQ0 = e -iQ0j._ ( 67 ) 
Then we derive the equation for determing the residue symmetry 

Oii T a)ijJj ®h + 2 Ji® Quh = 0. (68) 
Finally we obtain following conditions 



I = 0, a/ 3, 8 

2 
1 



?? = -02 



(69) 



We can rewrite the group element as 
where the new U(l) charge is 



(70) 



Q = i®Q-(t 3 + -Lt 8 ^) ® 1 (71) 

The Q charge in the CFL phase for u, d and s quarks with different colors are 
listed in Tab. J3)). One can check that the condensate made of (u\d 2 , u 2 d\), 
(d 2 s 3 , d 3 s 2 ) and (s 3 ui, s\u 3 ) is neutral because 

Qcond = Q(ui) + Q(d 2 ) = Q(u 2 ) + 

= Q(d 2 ) + Q(s 3 ) = Q(d 3 ) + Q(s 2 ) 

= Q(s 3 ) + Q( Ul ) = Q( Sl ) + Q(u 3 ) = 0. (72) 
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5.3 Symmetry pattern: spin-1 phases 

The condensate of spin-1 CSC transforms as an anti-triplet under SU(3) C 
transformation, and as a vector under spatial 50(3) rotation. It carries 
a Ub(l) charge of the baryon number. The symmetry group of the theory 
is then G = SU(3) C x £4(1) x 50(3). The condensate can be written as 
Ai = A ai J a (g) 7*, where 7* with i = 1,2,3 are Dirac matrices. The order 
parameter A a j is a complex matrix of 3 x 3 dimension and transforms as, 
A -► giAg 2 , where g x G 17(3) = SU(3) C x U b (l) and g 2 G 50(3). This can 
be seen by the transformation of the condensate, 

A ai J a (8> 7* A ai giJ a gl ® 327* 5iT = A ai J b (gi) ba ® {g^ijj 3 

= [(9l)baAai(g2)ij]Jb ® 7 J '- (73) 

Here 52 = exp^ie^^^^) in g 2 J % g 2 ' 1S the 50(3) representation in Dirac 
space and 52 = exp(i#fcjfc) in (g 2 )ijj J is the 50(3) representation in vec- 
tor space. By a suitable transformation, A can be cast into the following 
form, A = 5 + iA, where 5 and A are real symmetric and anti-symmetric 
matrices respectively. Since orthogonal rotation does not change the sym- 
mety of 5 and A, they can be after transformation written in the form 
5 = diag(Ai, A 2 , A3) and Ay = e^Wfc, where Aj,Wj are real numbers. The 
condensate then has the form, 



(74) 



Generally a complex matrix has 18 real parameters, 12 of which can be 
fixed by a transformation of symmetry group G. We can classify the order 
parameter A by its invariant subgroups of G, which can be found by looking 
for gi and g 2 satisfying gig^Ag 2 = A, which is equivalent to g\A = A and 
g\ Ag 2 = A. Since the matrix A defines a vector u7 = Wjej, A is invariant 
under a rotation g 2 = e lJt ^ G 50(3) with the axis in the direction of ~uj, 
i.e. 6 2 || To, i.e. g 2 Ag 2 = A. One can verify that the same rotation g 2 also 
makes 5 invariant if 5 is in the form, 

Sij - j^jij » ( 75 ) 

To see that we note that MiOJiUjjUJk = 0. The U(3) invariant subgroup whose 
element g± satisfies g±A = A is determined by the number of zero modes 





( A, 


iu; 3 


—iuj 2 


H 


-10)3 


A 2 


iuj\ 




V iu 2 




A3 
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order parameter 



Ai +ia 
-ia Ai 
A 2 , 












unbroken symmetry 



SO(2) v 



SO{2) v x U(1) L 



SO(2) v x U(1) L 



SU(2) L x SO(2) v x C/(1) L 



SO(3) v 



SU(2) L x SO(2) fl x U(1) L 



U(X)l 



SU(2) L x C/(1) L 



SO(2) v 



SO(2) v x C/(l)j 



name 



oblate 



cylindrical 



A 



CSL 



polar 



N 2 



axial 



planar 



Table 4: Classification and the symmetry breaking patterns of the order 
parameters, see Fig.l of [34, 36J. L/R denote the transformation matrix 
acting on the left/right hand side of the condensate, while SO(2) v denotes the 
transformation that acts on both right and left hand side of the condensate. 
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of A. The subgroup is then U(n) for the case of n zero modes of A. The 
classification of the order parameters is listed in Tab. [H see [34l 136] . 

We take a few examples. For the oblate phase, one can verify that the 
order parameter has 50(2) residue symmetry, g\Ag 2 = A, for g2(oblate) = 
e lJ 3 e 2 . Note that det A 7^ 0, so there is no residue symmetry inherited from 
7/(3). For the cylindrical phase, the residue symmetry is 50(2) x 7/(1), 
whose transformation is 

0i = diag(l,l, e "), g 2 = e u * e l (76) 

Note that there is one zero eigenvalue for A so the invariant subgroup is U(l). 
For the CSL phase, the residue symmetry is 5*0(3), whose transformation 
matrix obeying gj Ag 2 = A is any rotational matrix since A is a unit matrix, 
i.e. g 2 (CSL) = e iJ ^ with i = 1,2,3. 

A more convenient way of looking for invariant subgroup is by using 
infinitesimal transformation, 

A = g x Ag 2 = (l+i0 + iT a 0?)A(l + iJi0i) 

« A + id A + i9iT a A + %Q\AJi, (77) 

which requires 

8A + 0"T a A + 9\AJi = (78) 

for any subset of non- vanishing angles 0, Of and 8\. 

We analyze a few phases. The forms of order parameters A are listed 
in Tab [H (1) The oblate phase. Since J3A = AJ3, we obtain the only 
non-vansihing angles are 8\ and &\ and obey Of = — 2Q\. So we have verified 
that the invariant subgroup is SO(2), whose generator is 2T 2 (g) 1 — 1 (g> jj. 
(2) The cylindrical phase. We have 9f = -2^38 and 0? = -20f, which 
corresponds to residue symmetry SO(2) x U(l). The unbroken generators 
are 2T 2 ® 1 - 1 <g> J 3 T for 50(2) and + T 8 for 17(1). (3) The e phase. 

We obtain B\ = -26\ for 50(2) and {8 = ±0?,0f = ^0?) for 7/(1). The 
invariant generators are 2T 2 ® 1 - 1 O jj for 50(2) and | + T 2 + -^T 8 for 
77(1) (to see this one can verify that the generator behaves like a projector). 
(4) The A phase. It is a special case of the cylindrical phase. We obtain 
8\ = -28\ for 50(2) and 8\ = -2^38 for U{1), whose generators are 
2T 2 ® 1 — 1 (g) jj and — + T 8 respectively. Additionally we have d\ = 

-idf, 8\ = -±8f and 8\ = i{8\ - 8\) - 8\ for 577(2), whose generators are 
-1=( J=T 2 + T 8 ), ^(T 5 + T 6 ), and -k- (T 7 - T 4 ). It is more convenient to 
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transform the order parameter A which is Hermitian in Tab [4] by a unitary 
matrix hi into, 



It is easy to check that the unbroken generators are, T % with i = 1,2,3 
for SU{2), -i= + T 8 = &P l2 for 17(1), and 1 + J 3 for SO{2). Here we 
define P\i = diag(l,l,0) and P3 = diag(0,0,l). (5) The polar phase. The 
unbroken generators are T l with i = 1,2,3 for SU(2), + T 8 = ^Pu 
for U(l), and J3 for SO(2). For more comprehensive analysis of symmetry 
breaking features of the spin-1 pairings, see [34l 1361137] . 

5.4 More physics related to symmetry patterns 

There are a lot of interesting phenomena arising from symmetry patterns in 
color superconductivity. For example, the new U(l) charge corresponds to 
the new photon in some superconducting phases, which resembles and 
Z° bosons in electroweak theory, see Tab. ([5]) for the comparison between 
the electroweak theory and color superconductivity. One can do electroweak 
physics with the new photon [52j [53] . In presence of the new photon, mag- 
netic fields show a special behavior in the neutron star core [54]. One can 
study the transimission and reflection of the new photon at the interface be- 
tween different phases to reveal its connection to confinement [55]. Especially 
the new photon in the CFL phase can propagate in the color superconductor. 
Therefore the external magnetic field can penetrate the color superconduct- 
ing quark core modifying the gap structure [561 157] . and producing different 
phases with different symmetries and low-energy physics [58, 59]. Another 
interesting thing related to the residual SU(2) C symmetry in the 2SC phase 
is that one can construct SU(2) C effective theory, where the elementary exci- 
tation is the glueball [60J. Through decaying to photons these glueballs can 
be possible source of Gamma Ray Burst |61j. 

6 Deriving the gap equation 

In the above symmetry analysis it was demonstrated that a multitude of 
phases are possible when quarks form Cooper pairs. To decide which of 
those will be favorable at given T and one has to determine the phase 
with the largest pressure. Generally, the pressure will grow with the gain of 



A 




(79) 
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Table 5: Comparison to the standard model. 





Weinberg-Salam 


CSC with 77/ = 2,3 


group 


517(2) x 17(1) 
isospin, hypercharge 


SU(3) C x U(l) em 
color, electromagnetism 


field 


w 1 ,w 2 ,w 3 ,w 


Ai,--- ,A 8 ,A 
gluons, photon 


coupling 
constant 


9,sf 


9,e 


symmetry 
breaking 


SU{2) x 17(1) 
- U(l) em 


SU(3) C x U{l) em 


new field 


W + ,W~,Z°,A 


A lr -- ,A 7 ,A S ,A 


new coupling 
constant 


e = g' cos Ow 


e = e cos 


massless 
field 


A 


A 



condensation energy of the Cooper pairs. This in turn will depend on the 
magnitude of the respective color superconducting order or gap parameters. 
In this section we will derive the QCD gap equation and solve it in the weak 
coupling limit. 

We will use units in which 7j = c = &b = 1. We denote a 4- vector by 
capital letters, = (a?o,x), with x being a 3-vector of modulus |x| = x 
and direction x = x/x. For the summation over Lorentz indices, we use a 
notation familiar from Minkowski space, with metric g^ v = diag(+, —,—,—), 
although we exclusively work in compact Euclidean space-time with volume 
V/T, where V is the 3- volume and T the temperature of the system. Space- 
time integrals are denoted as J^ T dr f v d 3 x = J x . 

The partition function for QCD in absence of external sources reads 

Z = JvAexp{S A [A}} Z q [A\. (80) 

Here the (gauge-fixed) gluon action is 

1 



Sa[A] = j 



+ S gf [A] + S ghost [A], (81) 



where F£ v = d^A® — d v A a ^ + gf abc A b A^, is the gluon field strength tensor, 
S g f is the gauge-fixing part, and S g h os t the ghost part of the action. 
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The partition function for quarks in presence of gluon fields is 

Z q [A] = Jv$Vil>exp{S q [A,lj;,il>]} , (82) 

where the quark action is 

S q [A,$, V] = / m) {i^D x + ^70 - m) ^{X) , 

with the covariant derivative D x = d x — igAa{X)T a . Note that we focus on 
a single chemical poteintial and a single mass for simplicity. In fermionic sys- 
tems at nonzero density, it is advantageous to additionally introduce charge- 
conjugate fermionic degrees of freedom, 

il> c {X) = Ci?{X) 

^ c (x) = 4> T (x)c 

MX) = of c {X) 

MX) = ^(X)C (83) 

where the charge-conjugation matrix satisfies C~ l = C' = C T = —C, 
C~ l ^^C = —7^. We may then rewrite the quark action in the form 

S q [A,9,V] = \ f ^{X)g,\X,Y)^(Y) 

1 JX,Y 

+| f V{X)r£AXX)9(X) , (84) 
where we defined the Nambu-Gor'kov quark spinors 

and the free inverse quark propagator in the Nambu-Gor'kov basis 

C-i(XY) ( ^ G or\X,Y) \ 

y ° lA ' yj " V [G^-^Y) J ' 

with the free inverse propagator for quarks and charge-conjugate quarks 

[G±]-\X,Y) = (i 1 ,d x ± mQ -m)5^(X-Y) , 

The quark- gluon vertex in the Nambu-Gor'kov basis is defined as 

f t - ( 7 7° - 7 °tj ) • i 85 ) 
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The factors 1/2 in Eq. (l84|) compensate the doubling of quark degrees of 
freedom in the Nambu-Gor'kov basis. 

As we shall work in momentum space, we Fourier-transform all fields, 
as well as the free inverse quark propagator. Since space-time is compact, 
energy-momentum space is discretized, with sums (T/V) = T XmC^AO Sk- 
For a large volume V, the sum over 3-momenta can be approximated by an 
integral, (1/F)£ k ^ Jd 3 k/(2vr) 3 . For bosons, the sum over n runs over 
bosonic Matsubara frequencies = 2nirT, while for fermions, it runs over 
fermionic Matsubara frequencies = (2n + l)7rT. In our Minkowski-like no- 
tation for four- vectors, xq = t = —it, ko = —iu^J 1 . The 4-dimensional delta- 
function is conveniently defined as 5^{X) = 5(t)5^(x) = —i5(x°)8^ 3 \x). 
The Fourier transformations for the fields and free inverse quark propagator 
are given by 



vv K 




K,Q 



A£(X) = —=\] e - ip - x AZ(P). (86) 



The normalization factors are chosen such that fields in momentum space are 
dimensionless. The free inverse quark propagator is diagonal in momentum 
space, too, 

G- l (KO) - 1 ( ^oTW \ (4) 

where [Gq = 1^ ± // 7o - m. 

Due to the relations (|83|) . the charge-conjugate quark field in momentum 
space is related to the original field via tpc(K) = Ctp (~K) and tp c (K) = 
il) T {-K)C. The measure of the functional integral over quark fields can then 
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Figure 5: T2 is the sum of all QCD two-particle irreducible vacuum diagrams. 
Here, only the two-loop contributions are shown explicitly. Wavy lines cor- 
respond to gluon propagation and straight lines to quark propagation. 




+ 




+ 



be rewritten in the form 



A 

= M Yl dip(K)dTp(K)(hjj(-K)diJj(-K) 

(K,-K) 

= M' ]"[ d$(K)d^(K)dib c (K)diPc{K) 

(K-K) 

= M" Yl d9(K)d9(K) = V^Vy , (88) 

{K-K) 

with the constant normalization factors M, ,M' and M" . The last identity 
has to be considered as a definition for the expression on the right-hand side. 

Inserting Eqs. ([86]) - ([88]) into Eq. (|82j) . the partition function for quarks 
becomes 



1 



exp 



Here, we employ a compact matrix notation, 



= l^(K)[g^(K,Q)+gA(K,Q)]^>(Q) , (89) 



K,Q 



with the definition 



A(K,Q) = -j=t*A%{K-Q). 



(90) 



Since the binding of Cooper pairs is a non-perturbative effect we have to 
resort to a non-perturbative, self-consistent, many-body resummation tech- 
niques to calculate the gap parameter. For this purpose, it is convenient to 
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employ the CJT formalism [62]. The first step is to add source terms to the 
QCD action, 

S[A,V,V] = S A [A] + S q [A t %V] 

-> S[A,y,y] + JA + ^AKA 

+- (VH + + y]C$!) , (91) 

where we employed the compact matrix notation defined in Eq. (|89j) . J, H, 
and H are local source terms for the soft gluon and relevant quark fields, 
respectively, while K and K, are bilocal source terms. The bilocal source K, 
for quarks is also a matrix in Nambu-Gor'kov space. Its diagonal components 
are source terms which couple quarks to antiquarks, while its off-diagonal 
components couple quarks to quarks. The latter have to be introduced for 
systems which can become superconducting, i.e. where the ground state 
has a non-vanishing diquark expectation value, (^COip) / 0. One then 
performs a Legendre transformation with respect to all sources and arrives 
at the CJT effective action [621 E3] 

r[A,%%A,g] = S[A,V,V] -il^lnA^-iT^A^A-l) 

+ ^Tr q lng- l + ^Tr q (g^Q - l) 
+r 2 [4*,*,A,0] . (92) 

The tree-level action S[A, \P] now depends on the expectation values A = 
(A), ^ = (vp), and ^ = (^) for the one-point functions of gluon and quark 
fields. The traces denoted by Tr 9 run over gluonic degrees of freedom (i.e. 
over adjoint color) and momenta, while the traces Tr q run over quark degrees 
of freedom(i.e. fundamental color), flavor, Dirac and Nambu-Gor'kov indices 
as well as over momenta. The quantity 

! 6*S A [A] 
SA6A 

in Eq. (|92l) denotes the free inverse gluon propagator. A and g are the 
expectation values for the two-point functions, i.e., the full propagators, of 
gluons and quarks, respectively. The functional T2 is the sum of all two- 
particle irreducible (2PI) diagrams. These diagrams are vacuum diagrams, 
so they have no external legs. They are constructed from the vertices defined 
by S[A,ty,iSf], linked by full propagators A, g, cf. Fig. (0. 
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The expectation values for the one- and two-point functions of the theory 
are determined from the stationarity conditions 







5T _ 5F _ 5T 

W~6A~6g- (93) 

The first condition yields the Yang-Mills equation for the expectation value 
A of the gluon field. The second and third condition correspond to the Dirac 
equation for \& and respectively. Since and ^ are Grassmann- valued 
fields, their expectation values must vanish identically, \& = ($) = \P = 
(&) = 0. On the other hand, for the Yang-Mills equation, the solution A is 
in general non-zero but, at least for the two-flavor color superconductor con- 
sidered here, it was shown [64, 65J to be parametrically small, A ~ 4> 2 / \g 2 fi) , 
where (p is the color-superconducting gap parameter. Therefore, to sublead- 
ing order in the gap equation it can be neglected. 

The fourth and fifth condition (|93l) are Dyson-Schwinger equations for 
the gluon and quark propagator, respectively, 

A" 1 = Ao X + n, (94) 



7-1 



Go' + s . (95) 



where 



n ee -2^, (96) 

E = 20 (97) 

are the gluon and quark self-energies, respectively. The Dyson-Schwinger 
equation for the quark propagator (|95l) is a 2 x 2 matrix equation in Nambu- 
Gor'kov space, since 

T-t-l-l 



Gq1 - \ [G c -i- 



s - ( £ *- ) < w 

where S + is the regular self-energy for quarks and S~ the corresponding 
one for charge-conjugate quarks. The off-diagonal self-energies $ , the so- 
called gap matrices, connect regular with charge-conjugate quark degrees 
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Figure 6: The quark self-energy £ written in the Nambu-Gor'kov basis, 
cf. Eq. (|98l) . Single straight lines directing to the left /right correspond 
to ([G^] _1 + S T ) . Double straight lines to the left/right correspond to 
. Full/empty dots symbolizes the Cooper pair condensate ^ absorb- 
ing/emitting a quasiparticle pair. 




of freedom. A non-zero <J> corresponds to the diquark condensate. Only 
two of the four components of this matrix equation are independent, say 
+ S + and <J> + . Charge conjugation invariance of the action requires 
[G Y 1 + £~ = C{[G+] _1 + £+} T C _1 . Furthermore, the action has to be 
real valued, yielding the relation <I>~ = 7o[ < ^ + ]^7o- The quark propagator Q 
can be formally inverted, 

9 = ( 1+ % ) ' (99) 

where 

q ± = {[Gtr l +^ 

-<t>T ([g^]- 1 + ^y 1 ^ ± y 1 (ioo) 

is the propagator describing normal propagation of quasiparticles and their 
charge-conjugate counterparts in presence of diquark condensates, while 

s± = -([Gjj^ + ST)- 1 ^^ (101) 

describes the anomalous propagation of quasiparticles. It can be interpreted 
as the absorption (+) and the emission (— ) of a quasiparticle-pair by the 
condensate, for details, see Ref. [5]. 

In order to proceed one has to make an approximation for T2. We will 
restrict to two-loop 2PIs. Furthermore, pure gluonic loops are proportional 
to T 2 [66] and can be safely neglected in the limit of small temperatures. 
Consequently, we will consider only the first diagramm in Fig. [5] It is 
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an advantage of the CJT formalism, that truncating T2 will not destroy 
the self-consistency of the solution. The gluon self-energy is computed as 
II = — 2 5r2/(5A T , i.e., by cutting a gluon line in the first diagram of Fig. 
J!]). Thus, in our approximation II is simply a quark loop with Nambu- 
Gor'kov propagators Q. Accordingly, the quark self-energy £ = 2<5r 2 /5{7 T , 
is obtained by cutting a quark-line in the same diagram. Its Nambu-Gor'kov 
structure is shown in Fig. [6l The diagonal components T,^ 1 correspond to 
the ordinary self-energies for particles and charge-conjugate particles. The 
self-energies ^ symbolize the Cooper pair condensate connecting particles 
and charge-conjugate particles. In the following, also the term gap matrix 
will be used for $ + . Explicitly, the Nambu-Gor'kov components of £ read 

E+(K) = -g^Y j rT a Q + (Q)l u T b ^f v {K-Q) , (102) 
Q 

= -g 2 ^rTlG-(Q)7 u T^A^(K-Q) , (103) 
Q 

= 5 2 ^7^J~ + (Q) 7 ^A^(K-Q) , (104) 
Q 

®~{K) = 5 2 ^7^T a H-(Q) 7 % T A^(^-Q) • (105) 
Q 

Inserting these self-energies (and the corresponding one for the gluons) into 
the Dyson-Schwinger equations (|94l) and (|95l) one obtains a coupled set of 
integral equations which has to be solved self-consistently. In particular, 
the Dyson-Schwinger equations for the off-diagonal components ^ of the 
inverse propagator , i.e., Eqs. fll04f) and (|105l) . are the gap equations for 
the color-superconducting condensate. 

Even in the mean-field approximation applied here a completely self- 
consistent solution of the coupled set of Dyson-Schwinger equations is not 
feasible. Therefore, we will refrain from determining the gluon propagator 
self-consistently but resort to its hard-dense loop approximation [66]. Thus 
we will obtain a well-controlled, approximate solution of the gap-equation as 
it will be demonstrated in next section. 
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7 Quasi-particle excitations in color superconduc- 
tor 



In this section we analyse how the presence of a Cooper pair condensate 
affects the excitation spectrum of quasiparticles. To this end, the poles of the 
propagator Q have to be determined. The excitation spectrum depends on 
the specific symmetry broken pattern of the superconducing phase considered 
here, cf. Sec. [3] and 03 To this end, we first analyse the Dirac, color and 
flavor structure of the varios terms contained in Q, cf. Eq. (|100p . As 
we work in the ultra-relativistic limit, m = 0, we may expand the free 
(charge conjugate) quark propagator Gq in terms of positive and negative 
energy states using the projection operator A£ = (1 + e7o7 • k)/2, where 
e = + projects on (charge conjugate) quarks and e = — projects on (charge 
conjugate) antiquarks. One obtains for the respective free propagators 

[Gfrr^K) = 7 o£ [ko - (ek T (i)] K . (106) 

e=± 

Furthermore, [Gq ] _1 is diagonal in color and flavor space. To simplify the 
further analysis we anticipate that in weak coupling limit adopted here and 
to subleading order, the regular self-energy may be approximated by its 
one-loop approximation (67JEBIE9] in the leading logarithmic order (cf. Ref. 
[70] for result beyond the leading logarithmic order). In this approximation, 
the self-consistent propagators in Eqs. (|102|103p are replaced by free propa- 
gators. In particular, this amounts to neglecting the Cooper pair condensate 
in £ , so that is diagonal in color and flavor space. In Dirac space one 
finds [671 ES ES] 

E+(K) = £-(#) « 7 5 2 ^feoln-^- + i7r|fc |J , (107) 

where g = £?/(3\/2vr) and M 2 = (3vr/4)m2 (with m 2 g = JV/^). To sub- 
leading order it is sufficient to include only the real part of H^(K) [68]. Its 
effect is simply a shift of the quasiparticle poles, ko — > ko/Z(ko), where we 
introduced the quark wave function renormalization function 

/ M 2 \~ 1 
Z(k )^[l + g 2 \n^j . (108) 

The imaginary part of T^^K) can be shown to contribute beyond sublead- 
ing order [7l] and therefore will be neglected in the following. Physically, 
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Table 6: The structure of the matrices M.^ in various color-superconducting 
phases. The last two columns show the two eigenvalues A r of and their 
degeneracy d r (counting color-flavor degrees of freedom in the 2SC and CFL 
phases, and color-Dirac degrees of freedom in the three spin-one phases). 



phase 




Ai (di) 


A 2 (da) 


2SC 


75 0"2 J3 


1 (4) 


0(2) 


CFL 


75 I J 


4(8) 


1 (1) 


CSL 


J 


k + 7± (k) 




4(4) 


1 (8) 


polar 


J 3 


P + 7l(k) 




1 (8) 


0(4) 


A 


J 3 k x + ik v + j'l (k) + ij y ± (k) 


2(4) 


0(8) 



the imaginary part of E (K) gives rise to a finite life-time of the quasi- 
particles off the Fermi surface, which reduces the magnitude of the gap at 
sub-subleading order |71j . 

Turning to the structure of $ ± we introduce a compact notation following 
[30], 

$ + (K) = Y J <f> e (K)M ic Ai, (109) 

e=± 

where 4> e , the so-called gap function, is a scalar function of 4-momentum K. 
In general, the quantity VWk is a matrix in color, flavor, and Dirac space, 
which is determined by the symmetries of the color-superconducting order 
parameter. It can be chosen such that 

[-Mk,A k ]=0. (110) 

Note that the above is only valid in massless case. For the spin-0 phases 2SC 
and CFL the strucure of Ai^ has been determined in Sec. [H In Tab. [6] the 
corresponding structures for three spin-1 phases, the CSL, the polar and the 
A-phase, are shown. Note that these spin-1 phases are actually the mixtures 
of longitudinal and transverse phases, and that the transverse phases have 
the largest gaps. 

With representation (|109p for <J? + one may rewrite the second term on 
the r.h.s. of Eq. (fTOOll 

Q-G-Q+IG-}- 1 = |<m)| 2 £ k A k e , 

e=± 
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where we introduced 



L k = -foMlMk-yo . 



fill) 



Since L k is hermitian, it has real eigenvalues, A r , and can be expanded in 
terms of a complete set of orthogonal projectors, "P£, 



L k = KV } 



[112) 



In the five phases considered here, there are only two distinct eigenvalues 
and therefore two distinct projectors. The eigenvalues are also listed in Tab. 
El and the projectors can be expressed in terms of L k via 



1,2 



L k -A 



2.1 



A 



1,2 



A2,l 



(113) 



Note, that [L k , A k ] = [V^ 2 , A k ] = 0. Since the projectors V r A k form a basis 
in color, flavor, and Dirac space, inverting the term in the curly brackets in 
Eq. (llOOj) is straightforward, yielding 



i 



[k /Z(k )f 



-k,r' 



2 ' 



where 



ei ) M)=[(ek-^) 2 + \ r \ 
Inserting Eqs. (I109I114II into Eq. ([PUTT) yields 



-(K) = -^ 7o ^ k 7o^A l 



,21 1/2 



TO 



[fco/Z(fc )] : 



"k,r ^ 



2 ' 



(114) 



(H5) 



(H6) 



Obviously, the poles of the propagators Q + and H + are located at &o = 
±Z(ko) e kr ((/> e ) and the same result can be obtained for the propagators Q~ 
andH - . Eq. (11151) is the relativistic analogue of the standard BCS dispersion 
relation. The presence of the Cooper pair condenstate has generated a gap 
in the excitation spectrum of quasiquarks, e = +. Hence, even on the Fermi 
surface, k = fx, they carry a non-zero amount of energy, which is of the 
order of the gap function, e+ r = \/Ar </ ,+ - To excite (generate) a pair of 
a quasiparticle and a quasiparticle hole therefore costs at least the energy 
2\f\.(f) + . Such an excitation process can be interpreted as a break-up of a 
Cooper pair and the required energy as the Cooper pair binding energy. To 
investigate which quasiparticles exhibit a gap in their excitation spectrum 
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one has to determine the eigenvalues A r of the operator L^. For the phases 
considered here, has one or two non-zero eigenvalues. For A r = the 
respective quasi-particles remain ungapped. The degeneracy d r gives the 
number of excitation branches with the same gap, \/ cf)~^~ . 

Finally, for quasi anti-quarks, e = — , the dispersion relation hardly differs 
from the non-interacting case in the limit of large /i, k + (j,. 

In the next section, the gap equation (|104l) for the gap function (f) e {K) 
will be solved for quasi-particle excitations, e = +, at zero temperature in 
order to determine the magnitude of the gap in their excitation spectrum. 



8 Solving the gap equation 

To obtain the Dyson-Schwinger equation for the scalar gap function 



■(K) 



one inserts Eq. (I116P into Eq. (11041) . multiplies from the right with -M^A^, 
and traces over color, flavor, and Dirac space. The result is the so-called gap 
equation 



<P e '(Q) 



( 27r ) 3 fe[9o/^o)] 2 -[< s (0 e ')]- 



xU*{K- Q)T£f (k,q) 



where 



Tr 






Tr 







(117) 



fll8) 



The power counting scheme for the various contributions arising on the r.h.s. 
of the gap equation (|117|) is the same for all phases considered here. Due 
to the factor g 2 the integral must give terms proportional to (j)/g 2 in order 
to fulflill the equality. They contribute to the gap function at the largest 
possible order, O ((/)), and are therefore referred to as leading order contribu- 
tions. Terms of order <f>j g from the integral contribute at the order O(gcj)) to 
the gap function and are called subleading order corrections. Sub-subleading 
order corrections contributing at 0(g 2 4>) to the gap function have not been 
calculated yet and will be neglected in the following. 

In principle it is possible to proceed without deciding on one of the specific 
phases considered here. However, to keep the calculation as transparent as 
possible we restrict to the 2SC phase. The results for the other phases will 
be summarized at the end of this section. 



35 



As already mentioned in the previous section, the Dyson-Schwinger equa- 
tion for the gluon propagator will not be solved self-consistently. Instead we 
will adopt the gluon propagator in its HDL approximation. In pure Coulomb 
gauge they can be decomposed into a electric (or longitudinal) and a mag- 
netic (or transversal) part as [51 EE] 

A 00 (P) = A e (P) 
A 0l (P) = 

Atf(P) = A t (P) (<% - piPj) . 

with P = K — Q and the longitudinal and transverse propagators Ag ;t . Since 
the HDL propagators are diagonal in adjoint color space, A^ = 5 a bA^ v , the 
color indices can safely be suppressed. The explicit forms of the propagators 
A^t can be found in Ref. [51(66]. For our purposes the approximative forms 
of Agj given below are sufficient. As will be substantiated in the following 
the leading order contribution to the gap arises from almost static magnetic 
gluons with 

A t LDM (^)^ p6+ ^ 4a;2 e(M-p). (119) 

The term M 4 u> 2 arises from Landau-damping which dynamically screens 
magnetic gluons at the scale p ~ M 2//3 c<// 3 . Consequently almosr static 
magnetic gluons dominate the interaction between quarks at large distances. 
Subleading order corrections come from non-static magnetic gluons with 

Af SM (p)^\e(p-M) , (120) 

and from static electric gluons with 

A F00 * ~ 2 A 2 < ( 121 ) 

p z + irrig 

where the constant mass term m 2 g = Nf^r provides Debye-screening at the 
scale p ~ gfi. As the gap equation (11171) requires a self-consistent solution, 
we may proceed by first giving the known result and confirming it in the 
following. At zero temperature the 2SC gap function assumes for on-shell 
quasiparticles (e = +) with momenta k « p. the value 



2W/ 0/ uexp(-|:) , (122) 
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where 



b' 



256vr 4 



exp 



N f g 2 

7T 2 + 4 



5/2 



(123) 



The term in the exponent of Eq. (|122p was first determined by Son [72] via 
a renormalization group analysis taking into account the long-range nature 
of almost static, Landau-damped magnetic gluons, cf. Eq. (I119|) . Before 
that, a non-zero magnetic mass (analogous to the Debye-mass for electric 

gluons) was assumed leading to the standard BCS exponent, exp (— ^r^ j , 

where A 2 is the the typical momentum of the exchanged gluon and <? 2 /3 the 
effective coupling in the anti-triplet channel. This finding was crucial since 
in the correct exponential (I122p the parametric suppression becomes much 
weaker for small values of g, i.e. for large values of \x. The factor b in front 
of the exponential is generated by the exchange of non-static magnetic (I120|) 
and static electric (|12ip gluons [28j [73j [74] . The prefactor b' is due to the 
real part of the quark self-energy as given by Eq. (11081) [671168] . 
From Tab. (J6j) we read off 

M k = 7502J3 
and hence with Eqs. fjlll j 1 13|) one finds the projectors to be 



vl 



j. 



3 ! 



1 - n , 

with the eigenvalues Ai = 1 (4-fold) and A2 = (2-fold). This yields 



(k,q) 



T/' ee '(k,q) 



1+ee'k-q) , 



3 — ee'k • q . 

(ek - e'q) 2 



P 



1 + ee'k • q 



(124) 



and T£ u ' ee (k, q) = 0. Hence, the gapless quasi-particles corresponding to 7-^ 
do not enter the gap equation. Performing the Matsubara sum in Eq. (|117|) 
one then obtains for the gapped branch with the projector V\ to subleading 
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order 



2 



x tanh -3- > / (inn I 



' "if" I 1— 

{-2A^(p),L+[2Af SM (p) 
+A^ M ( eq ' + 4,p) 

+A^ M (e^-el P )]4 m }. (125) 

To subleading order the coefficients are given by t/q = «q = 2/3 and 
r/2 = — 772 = — 1/6- (Other phases than 2SC have different coefficients r/ 2 ^ 
and possibly have an additional excitation branch. To subleading order, 
however, the coefficients ry^ (with — 1 < m < 2) become simple numbers of 
order 0(1) for all phases considered here. It follows that the values of the 
respective gap functions differ only at subleading order [30].) To obtain Eq. 
(|125p several sub-subleading contributions are dropped. For example, the 
integration over q can be confined to a narrow interval of length 25 around 
the Fermi surface, 5 ~ g/j,, since the gap function peaks strongly around the 
Fermi surface. Furthermore, use was made of the fact, that the gap function 
must be an even function in K, 4> e {K) = cp e (—K). This property follows 
from the antisymmetry of the fermionic quark wave- functions. One has with 
ipc(K) = Ci> T {-K), tpc(K) = ip T (-K)C and C = -C" 1 = -C T 



\ T J.Ti 

K K 



Y,MK)C- 1 l<S> + {-K)] T C^(K). 



K 

Hence, the gap matrix must fulfill 

C<S> + {K)C- 1 = [$ + (-K)] T . 

Since in the 2SC case we have C^A^C -1 = [7sAi k ] T and [J3T2] 71 = J3T2, it 
follows that (fi e (K) = 4> e (—K). Moreover, the gap function was assumed to 
be isotropic in momentum space, 4> e {k) = 4> e (k) [30]. Finally, all imaginary 
contributions arising from the energy dependence of the gluon propagators 
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and the gap function itself are neglected. It follows, that the energy depen- 
dence of the gap function occuring in the final result has not been determined 
self-consistently. In Ref. [73J it is argued that imaginary contributions are 
negligible to subleading order. 

Using the well-known solution Eq. (I122|) for the 2SC gap function it is 
now straightforward to confirm that the various gluon sectors occuring in 
Eq. (|125p indeed contribute at the claimed orders of magnitude. For the 
static electric and the non-static magnetic gluons the integral over q gives 
for e' = + and m = 1 



/x-5 4 V <t> + 



2 5 Vln( =- 



26 



9 



In 



exp(-l/g) 



9^ 



(126) 



where for simplicity the gap function is moved out of the integral and the 
fact was used that i ~ j/i » ^ ~ //exp(— 1/g). Note that this result 
does not depend on the specific choice of S, as long as S is of the order 
5 ~ g c {ii a (f) b )°-+b with some numbers a/0 and b, c. The integral over p in 
combination with the factor 1/k ~ p, is of order 0(1). It follows that static 
electric and non-static magnetic gluons actually contribute at subleading 
order to the gap function. The above integral and the logarithm that follows 
from it also appear in normal BCS theory [75 j . Therefore we will refer to it 
as the BCS-logarithm. For e' = — , one has e~ ~ q + \x and the contribution 
corresponding to Eq. (I126|) is of order g 2 (j)5//j, ~ g 3 (j), i.e. beyond subleading 
order. 

In the case of almost static, Landau-damped magnetic gluons the integra- 
tion over q can be estimated as in Eq. (|126p . Choosing the term with e' = + 
and m = and canceling q ~ [i in the integral with the factor 1/k ~ 1/p 
at the front one obtains the BCS-logarithm. In order to render the total 
contribution of leading order, the p— integral must yield an additional factor 
1/g. Using the notation uj± 
the term 



e+iej^ one finds as the dominant contribution 



J k—q 



dp 



p' 



p % + M 4 i4 



6 



In 



M 6 + M 4 t4 



(k 

c4 



) 6 + M 4 i4 

1 

9 ' 



(127) 
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The last two steps are justified if one chooses for the cut-off d in the q integral 
5 = (M 2 ^) 1 / 3 . Consequently, in combination with the BCS-logarithm A126H 
this indeed is a leading order contribution. 

After having identified the leading and subleading order contributions, 
one evaluates the p integral in Eq. (I125|) and finds to subleading order 
(supressing all the superscripts e = e' = +) 

S 



(ek,fc) = 9 2 / d(q - n)Z{e^) 
Jo 



xtanh^llnl^f^l . (128) 

How to solve Eq. (|128p to subleading order is demonstrated explicitly in |68j . 
In the following, only some important steps are presented. As first proposed 
by Son [72] one approximates 

- @(e q , s - e k>r )ln ( — 





+©(ek,r - e q , s ) In 



£k,r 



Rewriting Eq. (1128(1 in terms of the new variables [68 



k- n + e k J ' 

y - Slnf^L) (129) 
\q- fJ- + e q J 

and differentiating twice with respect to x one can transform the gap equa- 
tion into the form of Airy's differential equation. Its solution can be written 

as 

<j>{x) = (f>f C F{x) , (130) 

where <^o SC is given by Eq. (11221) and F(x) is a combination of different 
Airy- functions [68]. The situation simplifies considerably when one neglects 
the effect of the normal quark self-energy, i.e. the factor Z(e q ) in Eq. (|128|> . 
Besides b' = 1 in Eq. (pC2g]) , one then finds F(x) = sin x [721 [73] , as the 
differential equation for <p simplifies to that of an harmonic oszillator. For 
\k — fj,\ ~ (f> one obtains x = ir /2 + 0(g), while for \k — fj,\ ~ g/J. it is x ~ 0(g). 
Hence, for momenta exponentially close to the Fermi surface, \k — fi\ ~ 0, 
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the on-shell gap function is 0(efc, k) = <Pq SC [1 + 0(g 2 )}, i.e. constant to 
subleading order. Otherwise, for example for momenta \k — fi\ ~ gfj,, one has 
already 4>(ek, k) = g~(j)Q SC . Hence, the on-shell gap-function is sharply peaked 
around the Fermi surface. 

The gap for on-shell quasi-particles at the Fermi surface in the 2SC phase 
as given by Eq. (|122l) is valid for asymptotic densities where g(ji) <C 1. At 
chemical potentials of physical relevance, however, one has g(fi) ~ 0(1)- 
Then sub-subleading terms which are corrections of order 0(g) to the pref- 
actor of the gap function in principle have to be accounted for. So far, they 
have not been calculated. Keeping this caveat in mind, it is still instructive 
to extrapolate the present result down to physical quark-chemical potentials 
of several hundred MeV. As explained in more detail in Ref. [5] one finds 
c/)f c ~ 10 MeV. 

As already mentioned before, the solution of the gap equation in all 
phases considered here [301 [34] can be written in a general form. For all 
considered phases it is found that differences arises only in the prefactor in 
front of the exponential, i.e. are only of subleading order [30] [34] : 

0o = 2bb' 0f ie~ d e~" exp , (131) 

where d = and q = ln((Ai) ai (A 2 ) a2 ) 1/2 for all spin-zero phases. Here the 
coefficients ai,a 2 are positive constants obeying a± + a% = 1. In spin-one 
phases d ^ and q are given by 

1 (ni A 9i i In A gi i + n 2 A 9>2 In A ?]2 ) g 
q = -i 

2 (niA 3i i + n 2 Ag i2 )q 

where n r is the degeneracy factor for the excitation branch. Generally the 
eigenvalues A r depend on the direction of the quark momentum, so the angu- 
lar average () is taken for all directions q. For all spin-one phases considered 
in this paper, d is about 5, it means e~ d ~ 10~ 2 — 10~ 3 , which strongly 
reduce spin-one gaps relative to spin-zero-gaps. Tab. (|7]) give an overview 
of the gaps of the different phases in units of the 2SC gap, ^o SC - 

Tab. ([7]) also lists the critical temperatures T c for all considered phases 
[301 1341 1681 173] , Generally the ratio of the critical temperatures to the zero 
temperature gaps on the Fermi surface is given by 

— = — e~ ? ~0.57e- ? (132) 
4>0 vr 

where 7 ~ 0.577 is the Euler-Mascheroni constant. An important feature is 
that in such single-gapped phases as the 2SC, polar or A phase, the ratio is 
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Table 7: The value of the gap function at the Fermi surface, (f>o in units of 
its value in the 2SC phase, and the critical temperature, in units of its value 
expected from BCS theory, and in units of the critical temperature in the 
2SC phase, T C 2SC . 



phase 


d 




T c /T c BtJS 




2SC 





l 


1 


1 


CFL 





2 -l/3 


2 l/3 


1 


CSL 


5 


2-1/8 e -rf 


2 l/6 


e~ d 


polar 


5 


e -d 


1 


e -d 


A 


21/4 


e -d 


1 


e -d 



the same as in BCS theory (the so-called BCS ratio T c /4>q = 0.57) [73]. For 
phases with more than two different excitations, such as the CFL or CSL 
phase, one has e _s ^ 1, which violats the BCS ratio [29, 30J, cf. the third 
column in Tab. (0). Another situation where the BCS ratio is violated is 
that the gap is anisotropic in momentum direction. Expressing the respective 
critical temperatures T c in terms of T C 2SC the factor e _? in Eq. (11311) cancels 
the factor in Eq. (fT32l) yielding T C /T 2SC = e~ d , which is listed in the 
forth column of Tab. ([7|). The transition temperature of a superconducting 
phase is normally regarded as the temperature where the gap disappears 
in the gap equation. This is justified in the mean field approximation. If 
one goes beyond the mean field approximation by including the fluctuation 
of the diquark fields, even above the transition temperature the 'gap' (the 
so-called pseudogap) is not vanishing [76l[77]. The transition temperature 
can also be modified by the gluon fluctuation which is beyond the mean field 
approximation [78] . 

We know that in the gap equation there are gluon propagators. If one 
changes the gauge parameter of the gluon propagator, the question arises if 
the gap would be invariant. The answer is yes under the following conditions 
[64] : (a) the gap equation is put on the mass-shell, and (b) the poles of 
the gluon propagator do not overlap with those of the quark propagator. 
However in the mean field approximation the gap depends on the gauge 
parameter at the sub-subleading level if one imposes the on-shell condition 
[791 ED]. 

In intermediate coupling regime, the superconducting phases can be in- 
vestigated within a Dyson-Schwinger approach for the quark propagator in 
QCD [ST], E2] . It was found that at moderate chemical potentials the quasi- 
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particle pairing gaps are several times larger than the extrapolated weak- 
coupling results. 

9 Debye and Meissner masses 

In this section we will address another two important energy scales in su- 
perconductivity: Debye and Meissner masses. The Debye mass for gluons/ 
photons characterizes the screening of the chromoelectric/electric charge in 
the medium, while the Meissner mass leads to expelling the chromomag- 
netic/magnetic field out of the type-I color/normal superconductor. We 
consider the weak coupling limit where the gap is parametrically smaller 
than the Debye or Meissner mass, so the coherence length (inverse of the 
gap) is much larger than the penetration length (inverse of the Meissner 
mass), which means in the weak coupling limit, all color superconductors 
are of type-I. 

The Debye and Meissner masses for gluons or photons can be calculated 
through gluon or photon polarization tensors and taking the static limit 
(po = 0, p — > 0). The polarization tensors depend on the diquark condensate 
through the quark propagators in the quark loop. The self-energies with 
full momentum dependence and the spectral densities of longitudinal and 
transverse gluons at zero temperature in color-superconducting quark matter 
was studied in detail in [83], [84], [85] . ^ n some color superconducting phases 
the U(l) charge which can propagate in the medium is the rotated or new 
photon. It can be found by symmetry analysis shown in previous section or 
by diagonalizing polarization tensors for both the photon and gluons. We 
can treat the photon in the same footing as eight gluons by enlarging gluon's 
color indices to denote the photon by a = 0. Diagonalizing the polarization 
tenors with respect to new 'color' indices a = 0, 1 — 8, we can find Deybe 
and Meissner masses of the rotated photon and gluon which are eigenvalues 
of the resulting mass matrices. The detail calculation is carried out in Ref. 
[86] . The results are listed in Tab. PfTUl). 

We distinguish between the normal-conducting and the superconducting 
phase. In the normal-conducting phase, i.e. for temperatures larger than 
the critical temperature for the superconducting phase transition, T > T c , 
Meissner mass is vanishing , i.e., as expected, there is no Meissner effect in 
the normal-conducting state. However, there is electric screening for tem- 
peratures larger than T c . Here, the Debye mass solely depends on the num- 
ber of quark flavors and their electric charge. We find for T > T c (with 
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Table 8: Zero-temperature Debye masses. All masses are given in units of 
N f fi 2 /(6ir 2 ), where N f = 2 in the 2SC phase, N f = 3 in the CFL phase, 
and Nf = 1 in the polar and CSL phases. We use the abbreviations C = 
(21 - 8 In 2)/54, a = (3 + 4 In 2) /27, and = (6 - 4 In 2)/9. 







m D,a 7 = 
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1 


2 


3 


4 


5 


6 


7 


8 


1-7 


8 





2SC 







3 3 2 








2e^ 


CFL 


3C5 2 





-2</3(eg 


4Ce 2 


polar 





3572 


3 3 2 








18g 2 e 2 


CSL 


3/V 


iag 2 


3/V 


3f3g 2 


3a 5 2 


3/33' 


3ag 2 


3/V 








18g 2 e 2 



Table 9: Zero-temperature Meissner masses. All results are given in the 
same units as the Debye masses in Table The abbreviations of Table [8] 
are used. 
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1-7 
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2SC 





.9 2 /2 


5 2 /3 





e 5 /(3V3) 


e 2 /9 


CFL 
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-2£eg/V3 


4Ce 2 /3 


polar 
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4g 2 e 2 


CSL 


/V 
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/V 
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/V 
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3q 2 e 2 
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Table 10: Zero-temperature rotated Debye and Meissner masses in units of 
Nf/j, 2 /(6ir 2 ) and mixing angles for electric and magnetic gauge bosons. The 
constants (, a, and (3 are defined as in Tables [8] and [9j 
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m 2 D,a~ ( = 
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where q n is the electric charge for the quark flavor labled by n. Consequently, 
the 9x9 Debye mass matrix is already diagonal. Electric gluons and electric 
photons are screened. 

The masses in the superconducting phases are more interesting. The 
results for all phases are collected in Tab. (|8|) (Debye masses) and Tab. (JU) 
(Meissner masses). The physically relevant, or "rotated", masses are obtained 
after a diagonalization of the 9x9 mass matrices. We see from these tables 
that all off-diagonal gluon masses, a, b = 1 — 8, as well as all mixed masses 
for a = 1 — 7 vanish. Furthermore, in all cases where the mass matrix is not 
diagonal, we find m| 7 = m^m^. We collect the rotated masses and mixing 
angles for all phases in Tab. (TTOj) . 

Let us first discuss the spin-zero cases, 2SC and CFL. In the 2SC phase, 
due to a cancellation of the normal and anomalous parts, the Debye and 
Meissner masses for the gluons 1,2, and 3 vanish. Physically, this is easy to 
understand. Since the condensate picks one color direction, all quarks of the 
third color, say blue, remain unpaired. The first three gluons only interact 
with red and green quarks and thus acquire neither a Debye nor a Meissner 
mass. We recover the results of Ref. [87]. For the mixed and photon masses 
we find the remarkable result that the mixing angle for the Debye masses is 
different from that for the Meissner masses, ^ / 6m- The Meissner mass 
matrix is not diagonal. By a rotation with the angle 9m, given in Tab. (TTO]) . 
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we diagonalize this matrix and find a vanishing mass for the new photon. 
Consequently, there is no electromagnetic Meissner effect in this case. The 
Debye mass matrix, however, is diagonal. The off-diagonal elements m 2 D §7 
vanish, since the contribution of the ungapped modes cancels the one of the 
gapped modes. Consequently, the mixing angle is zero, Od = 0. Physically, 
this means that not only the color-electric eighth gluon but also the electric 
photon is screened. Had we considered only the gapped quarks, we would 
have found the same mixing angle as for the Meissner masses and a vanishing 
Debye mass for the new photon. This mixing angle is the same as predicted 
from simple group-theoretical arguments. The photon Debye mass in the 
superconducting 2SC phase differs from that of the normal phase, which, for 
qi = 2/3 and q 2 = -1/3 is m 2 Dfn = 5 N f e 2 fi 2 /(6ir 2 ). 

In the CFL phase, all eight gluon Debye and Meissner masses are equal. 
This reflects the symmetry of the condensate where there is no preferred 
color direction. The results in Tab. JB|) and (|9j) show that both Debye and 
Meissner mass matrices have nonzero off-diagonal elements, namely mg 7 = 
m 2 8 . Diagonalization yields a zero eigenvalue in both cases. This means that 
neither electric nor magnetic (rotated) photons are screened. Or, in other 
words, there is a charge with respect to which the Cooper pairs are neutral. 
Especially, there is no electromagnetic Meissner effect in the CFL phase, 
either. Note that the CFL phase is the only one considered in this paper in 
which electric photons are not screened. Unlike the 2SC phase, both electric 
and magnetic gauge fields are rotated with the same mixing angle 9d = #m- 
This angle is well-known [Ml l88l 189] . 

Let us now discuss the spin-one phases, i.e., the polar and CSL phases. 
For the sake of simplicity, all results in Tab. J8]), Jl]), and (fTU]) refer to a 
single quark system, Nf = 1, where the quarks carry the electric charge q. 
After discussing this most simple case, we will comment on the situation 
where Nf > 1 quark flavors separately form Cooper pairs. The results for 
the gluon masses show that, up to a factor Nf, there is no difference between 
the polar phase and the 2SC phase regarding screening of color fields. This 
was expected since also in the polar phase the blue quarks remain unpaired. 
Consequently, the gluons with adjoint color index a = 1, 2, 3 are not screened. 
Note that the spatial z-direction picked by the spin of the Cooper pairs has 
no effect on the screening masses. As in the 2SC phase, electric gluons do 
not mix with the photon. There is electromagnetic Debye screening, which, 
in this case, yields the same photon Debye mass as in the normal phase. 
The Meissner mass matrix is diagonalized by an orthogonal transformation 
defined by the mixing angle. 



46 



In the CSL phase, we find a special pattern of the gluon Debye and 
Meissner masses. In both cases, there is a difference between the gluons 
corresponding to the symmetric Gell-Mann matrices with a = 1,3,4,6,8 
and the ones corresponding to the antisymmetric matrices, a = 2, 5, 7. The 
reason for this is, of course, the residual symmetry group 50(3) c+ j that 
describes joint rotations in color and real space and which is generated by a 
combination of the generators of the spin group 50(3) j and the antisymmet- 
ric Gell-Mann matrices, T2, T5, and T7. The remarkable property of the CSL 
phase is that both Debye and Meissner mass matrices are diagonal. In the 
case of the Debye masses, the mixed entries of the matrix, mp^, are zero, 
indicating that pure symmetry reasons are responsible for this fact (remem- 
ber that, in the 2SC phase, the reason for the same fact was a cancellation of 
the terms originating from the gapped and ungapped excitation branches). 
There is a nonzero photon Debye mass which is identical to that of the polar 
and the normal phase which shows that electric photons are screened in the 
CSL phase. Moreover, and only in this phase, also magnetic photons are 
screened. This means that there is an electromagnetic Meissner effect. Con- 
sequently, there is no charge, neither electric charge, nor color charge, nor 
any combination of these charges, with respect to which the Cooper pairs 
are neutral. 

Finally, let us discuss the more complicated situation of a many-flavor 
system, Nf > 1, which is in a superconducting state with spin-one Cooper 
pairs. In both polar and CSL phases, this extension of the system modifies 
the results in Tab. (JSj) , j9]), and (fTUI) . We have to include several different 
electric quark charges, q±,..., qN f , and chemical potentials, /ii, . . . , /ijv>- In 
the CSL phase, these modifications will change the numerical values of all 
masses, but the qualitative conclusions, namely that there is no mixing and 
electric as well as magnetic screening, remain unchanged. In the case of the 
polar phase, however, a many-flavor system might change the conclusions 
concerning the Meissner masses. While in the one-flavor case, diagonalization 
of the Meissner mass matrix leads to a vanishing photon Meissner mass, this 
is no longer true in the general case with arbitrary Nf. There is only a 
zero eigenvalue if the determinant of the matrix vanishes, i.e., if m 2 M §7 = 
"i-M,88 "iM,77- Generalizing the results from Tab. JHJ), this condition can be 
written as 



Consequently, in general and for fixed charges q n , there is a hypersurface in 
the Nf- dimensional space spanned by the quark chemical potentials on which 




(133) 
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there is a vanishing eigenvalue of the Meissner mass matrix and thus no elec- 
tromagnetic Meissner effect. All remaining points in this space correspond 
to a situation where the (new) photon Meissner mass is nonzero (although, 
of course, there might be a mixing of the eighth gluon and the photon). Eq. 
(|133p is trivially fulfilled when the electric charges of all quarks are equal. 
Then we have no electromagnetic Meissner effect in the polar phase which 
is plausible since, regarding electromagnetism, this situation is similar to 
the one-flavor case. For specific values of the electric charges we find very 
simple conditions for the chemical potentials. In a two flavor system with 
qi = 2/3, (/2 = — l/3, Eq. (1133(1 reads /ufjU^ = 0. In a three-flavor system 
with qi = —1/3, q2 = —1/3, and (73 = 2/3, we have (/if + ^l)/- 4 ! = 0- Con- 
sequently, these systems always, i.e., for all combinations of the chemical 
potentials /%, exhibit the electromagnetic Meissner effect in the polar phase 
except when they reduce to the above discussed simpler cases (same electric 
charge of all quarks or a one-flavor system). 

In the weak coupling limit, the gap values are parametrically smaller than 
the Meissner masses, therefore the coherence length (inverse of the gap) is 
much larger than the penetration length (inverse of the Meissner mass). 
So the spin-1 color superconductor is a type-I superconductor in the weak 
coupling limit [32]. In the intermediate densities, this statement might not 
be valid. A natural way of studying superconductivity in external magnetic 
field is to make use of the Ginzburg-Landau effective theory, see e.g. Ref. 
[901 ED ESI E31 EH E51 [96]. An introduction to the Ginzburg-Landau theory 
can be found in Sec. [HI 

10 General effective theory for color superconduc- 
tivity 

As we mentioned before, quark matter at small temperature T and large 
quark chemical potential // is a color superconductor. In a color super- 
conductor, there are several energy scales: the quark chemical potential fj,, 
the inverse gluon screening length m g ~ g/x, and the superconducting gap 
parameter <\>. In weak coupling these three scales are naturally ordered, 
(ft <C gfJ- <C \x. The ordering of scales implies that the modes near the Fermi 
surface, which participate in the formation of Cooper pairs, can be con- 
sidered to be independent of the modes deep within the Fermi sea. This 
suggests that the most efficient way to compute properties such as the color- 
superconducting gap parameter is via an effective theory for quark modes 
near the Fermi surface. Such an effective theory has been originally pro- 
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posed by Hong [97] and was subsequently refined by others [971 1981 1991 HOOj . 

In order to construct such an effective theory, we pursue a different venue 
and introduce cut-offs in momentum space for quarks, A q , and gluons, A g . 
These cut-offs separate relevant from irrelevant quark modes and soft from 
hard gluon modes. We then explicitly integrate out irrelevant quark and 
hard gluon modes and derive a general effective action for hot and/or dense 
quark- gluon matter [1011 1102j . 

We show that the standard HTL and HDL effective actions are contained 
in our general effective action for a certain choice of the quark and gluon cut- 
offs A q , A g . We also show that the action of the high-density effective theory 
derived by Hong and others is a special case of our general effective action. 
In this case, relevant quark modes are located within a layer of width 2A q 
around the Fermi surface. 

The two cut-offs, A q and A g , introduced in our approach are in principle 
different, A q ^ A g . We show that in order to produce the correct result 
for the color-superconducting gap parameter to subleading order in weak 
coupling, we have to demand A, < j/i « A 9 < fi, so that A q /A g ~ g <C 1. 
Only in this case, the dominant contribution to the QCD gap equation arises 
from almost static magnetic gluon exchange, while subleading contributions 
are due to electric and non-static magnetic gluon exchange. 

The color-superconducting gap parameter is computed from a Dyson- 
Schwinger equation for the quark propagator. In general, this equation cor- 
responds to a self-consistent resummation of all one-particle irreducible (1PI) 
diagrams for the quark self-energy. A particularly convenient way to derive 
Dyson-Schwinger equations is via the Cornwall-Jackiw-Tomboulis (CJT) for- 
malism. In this formalism, one constructs the set of all two-particle irre- 
ducible (2PI) vacuum diagrams from the vertices of a given tree-level action. 
The functional derivative of this set with respect to the full propagator then 
defines the 1PI self-energy entering the Dyson-Schwinger equation. Since it 
is technically not feasible to include all possible diagrams, and thus to solve 
the Dyson-Schwinger equation exactly, one has to resort to a many-body 
approximation scheme, which takes into account only particular classes of 
diagrams. The advantage of the CJT formalism is that such an approxima- 
tion scheme is simply defined by a truncation of the set of 2PI diagrams. 
However, in principle there is no parameter which controls the accuracy of 
this truncation procedure. 

The standard QCD gap equation in mean-field approximation studied 
in Refs. [271 EH EH EH Q3EJ EH] follows from this approach by including 
just the sunset-type diagram which is constructed from two quark-gluon ver- 
tices of the QCD tree-level action. We also employ the CJT formalism to 
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derive the gap equation for the color-superconducting gap parameter. How- 
ever, we construct all diagrams of sunset topology from the vertices of the 
general effective action derived in this work. The resulting gap equation 
is equivalent to the gap equation in QCD, and the result for the gap pa- 
rameter to subleading order in weak coupling is identical to that in QCD, 
provided A g < gfi <C A g < /j,. The advantage of using the effective theory 
is that the appearance of the two scales A g and A g considerably facilitates 
the power counting of various contributions to the gap equation as compared 
to full QCD. We explicitly demonstrate this in the course of the calculation 
and suggest that, within this approach, it should be possible to identify the 
terms which contribute beyond subleading order to the gap equation. Of 
course, for a complete sub-subleading order result one cannot restrict oneself 
to the sunset diagram, but would have to investigate other 2PI diagrams as 
well. This shows that an a priori estimate of the relevance of different con- 
tributions on the level of the effective action does not appear to be feasible 
for quantities which have to be computed self-consistently. 



11 Ginzburg-Landau theory for CSC 

In this section we will introduce the Ginzburg-Landau theory for CSC. The 
starting point is the thermodynamic potential tt, which can be expanded 
in the second and fourth power of the diquark condensate in the vicinity of 
critical temperature. We have 

n = --^logdetS- 1 = -^Trlog^- 1 , (134) 

where Q~ l is the inverse of the full propagator and the trace runs over 
Nambu-Gorkov, momentum, Dirac, color and flavor space. We can expand 
O to the fourth power of E, 

n = -^THog^+E) 
= -llTHog^l + SoE)] 

-^TVlog^-^Trloga + aoE). (135) 
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Here the free part Q Q 1 is given by 
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where we only keep the positive energy excitation in the last line. Here 



the energy projectors are defined by = ^ 
Vk 2 + m 2 . The free propagator is then 
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The condensate part is 
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where $ = 7 o[ ( l )+ ]^7o- To simplify the last term in Eq. (|135p . we use 
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Eq. (1135(1 can be expanded as, 
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where the odd terms are vanishing due to zero trace in Nambu-Gorkov space. 

First we evaluate the quartic term of the condensate. Since we consider 
the small momentum expansion, we assume that £ is constant in the quartic 
term. Then we have 
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where we have used ko = i(2n + l)nT and 



[(2n + l) 2 vr 2 + x 2 ] 2 4x 3 



tanh 



x x/2 
2 ~ cosh(x/2) 



About the momentum integral in Eq. (|141|) . we see that the integrand except 
k 2 is centered at /i, if we limit the momentum range around fj, we can replace 
k 2 with fi 2 and pull it out of the integral. The first term inside the curly 
bracket dominates the momentum integral, then we have 
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The dependence on color superconducting phases contains in the average 
over momentum direction. 

Now let's evaluate the quadratic term, 
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where we have kept the momentum dependence of condensate. To the leading 
order the trace term invloves an angular integral over k. We can make Taylor 
expansion with respect to q for the last line of Eq. fll43|) . 
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Then the last line of Eq. (|143j) becomes 



kg 

2E k 



2£ fc ^ 8£ fe 



f kg 



+ 



+ 



(fc - (E k - /x)) 2 (fco + (E k - /i)) 2 fc 2 - (E k - jj) 



• (145) 



The second and third terms inside the square brackets are vanishing when 
carrying out integral over k since the term with is odd in k and the term 

2 

with is proportional to (E k — fi). Then the quadratic term becomes, 
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Note that the constant term inside the curly brackets gives the quadratic 

k 2 a 2 

term in the condensate. The term proportional to g, a gives the gradient 

k 

term ~ |V$| 2 . 

By inserting the concrete forms of the diquark condensates into Eq. (11421) 
and (11461) one can have the Ginzburg-Landau free energy from which one can 
study a variety of problems [90l [9TI [92l [93l f94l [95l [96]. 



12 BCS-BEC crossover in a boson-fermion model 

In weak couplings the Cooper pairing is well described within Bardeen- 
Cooper-Schrieffer (BCS) theory. Cooper pairs are typically of a size much 
larger than the mean interparticle distance. The picture changes for suffi- 
ciently large interaction strengths. In this case, Cooper pairs become bound 
states, and superfluidity is realized by a Bose-Einstein condensation (BEC) 
of molecular bosons composed of two fermions. A crossover between the 
weak-coupling BCS regime and the strong-coupling BEC regime is expected. 
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In order to describe the crossover from BCS to BEC in color supercon- 
ductivity, one can use Nambu-Jona-Lasinio model [1051 11061 11071 11081 11091 
1 11 01 1 1 111 11121 [TT3] . a purely fermionic model. Here we shall not consider a 
purely fermionic model which may describe this crossover as a function of 
the fermionic coupling strength. We rather set up a theory with bosonic and 
fermionic degrees of freedom. Here, fermions and bosons are coupled through 
a Yukawa interaction and required to be in chemical equilibrium, 2/i = 
where fi and /z& are the fermion and boson chemical potentials, respectively. 
We treat the (renormalized) boson mass m^ r and the boson-fermion coupling 
g as free parameters. Then, tuning the parameter x = — (m 2 r — ^)/(4g 2 ) 
drives the system from the BCS to the BEC regime. The fermionic chemical 
potential shall be self-consistently determined from the gap equation and 
charge conservation. This picture is inspired by the boson-fermion model 
of superconductivity considered in Ref. [114j . which has been used in the 
context of cold fermionic atoms. It also has possible applications for high- 
temperature superconductivity. For simplicity, we only discuss the evalua- 
tion of the model in a mean-field approximation (MFA) |115j . for the case 
beyond MFA, see Ref. [TTfi] . 

We use a model of fermions tp and composite bosons cp coupled to each 
other by a Yukawa interaction. With mean field approximation, the La- 
grangian can be written as 

c = hzs-^ + c^l-mlM 2 

+\(d t -i^\ 2 -\V<p\ 2 -m 2 b \tp\ 2 , (147) 
S-HP) = ( P ^ + ^- m *™<P V (148) 

Here charge conjugate spinors are defined by ific = Cip and ip c = ip T C with 

„2 „,0 __j tvt—u.. n— i : j„c_~ u..,t, _ ( ^ 



C = «7 7 , and Nambu-Gorkov spinors are defines by\P = I ^ j , ^ = 

{ipiipo), 5 _1 is the inverse fermion propagator which reads in momentum 
space, the fermion (boson) mass is denoted by m (mi,). We choose the bo- 
son chemical potential to be twice the fermion chemical potential, //& = 2/i. 
Therefore, the system is in chemical equilibrium with respect to the conver- 
sion of two fermions into one boson and vice versa. This allows us to model 
the transition from weakly-coupled Cooper pairs made of two fermions into 
a molecular difermionic bound state, described as a boson. The interac- 
tion term accounts for a local interaction between fermions and bosons with 
coupling constant g. In order to describe BEC of the bosons, we have to 
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separate the zero mode of the field tp. Moreover, we shall replace this zero- 
mode by its expectation value <f> = (</?o) and neglect the interaction between 
the fermions and the non — zero boson modes. This corresponds to the 
mean-field approximation. 

We can compute the partition function 



= J [d¥][dtf][dp][dp*]exp 



l/T 



drd 3 xC 



(149) 



to obtain the thermodynamic potential density = —T/V InZ, where T is 
the temperature and V is the volume of the system. After performing the 
path integral and the sum over Matsybara frequencies one may get 
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with quasi-particle energy for fermions (e = +1) and antifermions (e = 
— 1) denoted by e e k = \J (e^o — eu) 2 + A 2 , withe^o = Vfe 2 + m 2 , and the 

(anti)boson energy denoted by uif, = \ k 2 + 
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eu b . 



Then the charge conservation equation and gap equation can be derived 
by n = — and = which can be writen as n = np + wo + «b, with 
up, no, represent the fermionic, condensate bosonic and thermal bosonic 
contribution to the total partical density respectively, and for non-zero con- 
densation, the gap equation reads 
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with crossover parameter defined by x = 
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This parameter x can be varied from negative values with large modulus 
(BCS) to large positive values (BEC). In between, x = is the unitary limit. 
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Solving these two coupled equations to determine the gap A and the chemical 
potential fx as functions of the crossover parameter x and the temperature T 
at fixed effective Fermi momentum pp, fermion mass m, and boson- fermion 
coupling g (in numerical calculation, we fix ^ = 0.3, ^ = 0.2, g = 4, with 
A the cutoff in the momentum integrals). The solution [A(x, T), n(x ,T)] 
can then, in turn, be uesed to compute the densities of fermions and bosons 
in the x — T plane. 

The numerical results for the solution of the coupled equations are shown 
in Fig. [Jj The left panel shows the fermion chemical potential and the 
gap Ao as functions of x. In the weak-coupling regime (small x) we see that 
the chemical potential is given by the Fermi energy, = ep. For the given 
parameters, ep/A ~ 0.36. The chemical potential decreases with increasing 
x and approaches zero in the far BEC region. The gap is exponentially small 
in the weak-coupling region, as expected from BCS theory. It becomes of 
the order of the chemical potential around the unitary limit, x = 0, and 
further increases monotonically for positive x. In the unitary limit, we have 
fi/ep ~ 0.37, while in nonrelativistic fermionic models fi/ep ~ 0.4 — 0.5 was 
obtained [iTH filTH [TT9] . 

The corresponding fermion and boson densities are shown in the right 
panel of Fig. [Jj These two curves show the crossover: at small x all 
Cooper pairs are resonant states, which is characterized by a purely fermionic 
densityn = rip; at large x, on the other hand, Cooper pairs are bound states 
and hence there are no fermions in the system. The charge density is rather 
dominated by a bosonic condensate, n = tiq. The crossover region is located 
around x = 0. We can characterize this region quantitatively as follows. We 
write the boson mass as nib = 2m — -Ebind- Then, a bound state appears for 
positive values of the binding energy -E^ind) i.e., for 2m > mj. 

Let us discuss the quasi-fermion and quasi-antifermion excitation energies 
given by and . Inserting the numerical solutions for \x and A into these 
energies results in the curves shown in Fig. [8l These excitation energies show 
that, for large values of x, quasi-fermions and quasi-antifermions become 
degenerate due to the vanishingly small chemical potential. Because of the 
large energy gap, we expect neither quasi-fermions nor quasi-antifermions to 
be present in the system. 

The condensate and the fermion and boson density fractions at fixed 
crossover parameter x/xq = 0.2 are shown in Fig. [9j At the left end, T = 0, 
one recovers the results shown in Fig. [Jjat the particular value x/xq = 0.2, 
at the point T/T c = 1 the second-order phase transition manifests itself in a 
kink in the density fractions and a vanishing condensate. Below T c we ob- 
serve coexistence of condensed bound states, condensed resonant states, and, 
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Figure 7: Crossover at zero temperature from the BCS regime (small x) to 
the BEC regime (large x). Left panel: fermion chemical potential fio (blue 
dotted) and gap Ao (red dashed) in units of effective Fermi energy ep. Right 
panel: condensate fraction (red solid), fermion fraction (blue solid). 
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Figure 8: Fermion and antifermion excitation energies et and for three 

different values of the crossover parameter x/xq at T = as a function of 

the momentum k (both e e k and k are given in units of ef)- In the BCS 

regime (left panel) the energy gap is small and the fermion excitations are 

well separated from antifermion excitations. Both excitations approach each 

other in the unitary regime (middle panel), and become indistinguishable in 

the far BEC regime (right panel). Note in particular that the minimum of 

the antiparticle excitation is not a monotonic function of x. 
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Figure 9: Density fractions in the crossover regime at fixed x/xq = 0.2 
as functions of temperature: condensed bosons (red solid), fermions and 
uncondensed bosons (blue dotted and red dashed, respectively). 
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for sufficiently large temperatures, uncondensed bound states. We obtain 
thermal bosons even above the phase transition. They can be interpreted as 
"preformed" pairs, just as the uncondensed pairs below T c . This phenomenon 
is also called "pseudogap" in the literature [77]. It suggests that there is a 
temperature T*(x) which marks the onset of pair formation. This temper- 
ature is not necessarily identical to T c . In the BCS regime, T*(x) = T c (x), 
while for x > 0, T*(x) > T c (x). Of course, our model does not predict any 
quantitative value for T* because thermal bosons are present for all tem- 
peratures. Therefore, we expect the model to be valid only for a limited 
temperature range above T c . 

With the mean field appoximation, the boson-fermion model can be em- 
ployed to describe the BCS-BEC crossover. An important property of the 
model is the coexistence of weakly-coupled Cooper pairs with condensed 
and uncondensed bosonic bound states. In the crossover regime as well as 
in the BEC regime, strongly-bound molecular Cooper pairs exist below and 
above the critical temperature T c . Above T c , they are all uncondensed ("pre- 
formed" Cooper pairs) while below T c a certain fraction of them forms a 
Bose-Einstein condensate. In contrast, in the BCS regime, pairing and con- 
densation of fermionic degrees of freedom (in the absence of bosons) both 
set in at T c . 

It is straightforward to extend our boson-fermion model to two fermion 
species with cross-species pairing [1201 [39] . This allows us to introduce a 
mismatch in fermion numbers and chemical potentials which imposes a stress 
on the pairing |121U122] , This kind of stressed pairing takes place in a variety 
of real systems. For example, quark matter in a compact star is unlikely to 
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exhibit standard BCS pairing in the color-flavor locked (CFL) phase, i.e., 
pairing of quarks at a common Fermi surface. The cross-flavor (and cross- 
color) pairing pattern of the CFL phase rather suffers a mismatch in chemical 
potentials in the pairing sectors bu — rs and bd—gs (r, g, b meaning red, green, 
blue, and u,d,s meaning up, down, strange). This mismatch is induced by 
the explicit flavor symmetry breaking through the heaviness of the strange 
quark and by the conditions of color and electric neutrality. Our system 
shall only be an idealized and simplified model of this complicated scenario. 
However, as in the previous sections, we shall allow for arbitrary values of 
the crossover parameter and thus model the strong coupling regime of quark 
matter. We shall fix the overall charge and the difference in the two charges. 
This is comparable to the effect of neutrality conditions for matter inside a 
compact star, which also impose constraints on the various color and flavor 
densities. Our focus will be to find stable homogeneous superfluids in the 
crossover region and, by discarding the unstable solutions, identify parameter 
values where the crossover in fact becomes a phase transition. 

We denote the average chemical potential and the mismatch in chemical 
potentials by Ji = /J1 ^" /J2 , <5/x = Ml ~^ 2 . Then, the bosonic chemical potential 
is fj,b = 2JI . The thermodynamic potential differs from the one-fermion case 
in the dispersion relation for the fermions, 



n 



mi 



" V 

E 



A 2 + 



(2vr) 3 

+Tln [l + e-^-W]} 
With coupled gap equaiton (0 



{4 



+ Tln 



(153) 



an x 

dA' 
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and fixed sum (n = ni+ri2 
difference (5n = n\ — ri2 = — Jj^) °f the particle number density equations, 
the variables /I,5fi andA can be solved out. 

In order to check the gapless states for their stability, we have to compute 
the number susceptibility matrixxij = ^ , i = 1,2. Note that we fix n and 
5n (or equivalently rt\ and ni) in our solution. Hence \ can be regarded as 
measuring the response of the system to a small perturbation away from this 
solution. In particular, a stable solution requires the mismatch in density 
to increase for an increasing mismatch in chemical potentials. Therefore, a 
negative eigenvalue of this 2x2 matrix indicates the instability of a given 
solution. 
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Figure 10: The phase diagram in the plane of the crossover parameter x 
and the density difference Sn/n. Shaded areas have unstable homogeneous 
solutions with negative number susceptibility. NS denotes "normal state"; in 
this region, no solution for the gap equation is found. gSF denotes "gapless 
superfluid"; in this region a stable gapless superfluid state is found with two 
different Fermi surface topologies, divided by the right dashed-dotted line. 
Roman numbers II, III and IV denote three types of Fermi surface topologies 
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We finally present a phase diagram in Fig. [TU] for arbitrary (positive) 
values of 5n/n. Since we do not consider spatially inhomogeneous phases, 
this phase diagram is incomplete. Its main point is to identify regions where 
homogeneous gapless superfluids may exist. We find that for sufficiently 
large mismatches, 5n/n > 0.02 there is a region where no solution with 
nonzero A can be found. We see that the region of stable superfluids shrinks 
with increasing 5n/n, as expected. Note that the horizontal axis5n/n = is 
not continuously connected to the rest of the phase diagram. For vanishing 
mismatch in densities a stable, fully gapped superfluid exists for all x. One 
should thus not be misled by the instability for arbitrarily small mismatches. 

We conclude with emphasizing the two main qualitative differences to 
analogous phase diagrams in nonrelativistic systems: (i) within the stable 
region of homogeneous gapless superfluids there is a curve that separates two 
different Fermi surface topologies; this is the right dashed-dotted line in Fig. 
[TUl (ii) for large x the gapless superfluid becomes unstable even in the far 
BEC region; this is the shaded area on the right side in Fig. [TUl 
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13 Phenomenological implication: neutrino processes 
in compact stars 

Generally the color superconducting quark matter cannot be easily produced 
in laboratory. The most possible place for the color superconductor is in the 
core of some compact stars where the high baryon density is realized. In 
recent years, people have been looking for the connection of color supercon- 
ductivity to observations. The compact stars are born from supernova ex- 
plosions where neutrino processes play crucial roles. Looking for evidences of 
color superconductivity from neutrino processes is a natural choice. One pos- 
sible way is to study the compact star cooling through the neutrino emission, 
which is the major cooling process for proto-neutron stars and compact stars 
with relatively high internal temperature in the neutrino dominant stage. In 
this stage star bodies are transparent to neutrinos while opaque to photons. 
When the compact stars cool down by the neutrino emission, the photon 
emission gradually takes over for the later cooling. The compact star cool- 
ing models depend heavily on models for the neutron star interior especially 
the core. The star core could be nucleon matter, hyperon matter, pion or 
kaon condensate, or quark matter. There is huge amount of literature about 
compact star cooling since it was first studied by Tsuruta and Cameron in 
1966 [123]. For a recent review of the neutron star cooling, see, e.g. [124]. 
For a review of neutron star evolution, see e.g. [125J. A recent systematic 
work in neutron star cooling based on hadronic models can be found in Ref. 

Calculating neutrino prcocesses in normal quark matter is part of effort 
to describe the cooling behavior of compact stars with quark matter cores. 
Roughly speaking there are three kinds of the neutrino prcocesses in quark 
matter: the Urea, the modified Urea and the bremsstrahlung processes, see 
Fig. (TTTT) . The first systematic and complete calculation for the Urea process 
in normal quark matter was done by Iwamoto addressing the importance of 
Fermi liquid behavior for relaxing the strict constraint of energy-momentum 
conservation |127| . The recent work by Schafer and Schwenzer |128] extends 
Iwamoto's work by taking into account of the non-Fermi liquid behaviors 
[67 tl68l[T29lll30] . For the neutrino emission in fully gapped superconducting 
phases, the Urea process is exponentially suppressed |13 11 11321 11331 1134] . 
The possible relevant processes for the neutrino emission would be from the 
Goldstone mode |135] or from the pair breaking and recombination processes 
[13611137] . For gapless phases the Urea processes are not suppressed and like 
in the normal quark matter [j~38l [i~39l \TiO\ H4T] . 
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Figure 11: Neutrino processes in quark matter, (a) Urea; (b) Modified Urea; 
(c) Bremsstrahlung 
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Figure 12: Ordering along the close-time path in Green function. 
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In this section we will use close-time-path method to derive the Kadnoff- 
Baym equation. Based on the Kadnoff-Baym equation, we will obtain the 
kinetic equation for the neutrino emissivity for the normal and superconduct- 
ing quark matter. There are other applications of close-time-path method, 
see, e.g. Ref . [1421 H43j in the transport equation of the quark gluon plasma. 
For the derivation of Kadnoff-Baym equation and kinetic equation for other 
neutrino processes in compact star, see, e.g. Ref. [1441 1145] , If the neutron 
star has a color superconducting quark matter core which is in the so-called 
spin-1 A-phase, we will show that the neutrino emission is asymmetric in 
space. 



13.1 Covariant Kadnoff-Baym equation 

The two-point Green function for a fermion can be defined in the following 

way 

*G(l,2) = (T^ff(l)^(2)]), 

where 'V and '2' denote all indices including spatial coordinates of the two 
fields. T is the causal time-ordering operator which orders an operator at 
an earlier time to the right of an operator at a later time. The field opera- 
tors ipii(t) and i/jjj(t) in the Heisenberg picture are related to those in the 
Schodinger picture by 

1) H {t) = exp[iH(t-t )]^ s (t)exp{-iH{t-t )] 
tp H (t) = exp\iH{t - t )]ip s (t) exp[-iH(t - to)], 



where TL is the Hamiltonian of the system. 

In the same way, we can define the two-point Green function on the 
close-time path C as shown in Fig. (jl~2|) . 



*<?(!, 2) 



(7bhMl)^(2)]) 



exp 



c 



(154) 
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where £'j is the interacting part of the Lagrangian in the interaction picture 
and Tq is the ordering operator along the close-time path C . G(l, 2) can be 
written as 

G(l, 2) = e c (h - t 2 )G>(l, 2) + Q c (t 2 - ti)G<(l, 2), (155) 

where G >,< are correlation functions. We also need the retarded and ad- 
vanced Green functions as follows 

G R (1,2) = G(l,2)-G < (l,2) = O c (t 1 -t 2 )[G > (l,2)-G < (l,2)} 
G A (1,2) = G(l,2)-G > (l,2) = @ c (t 2 -t 1 )[G < (l,2)-G > (l,2)}. 

The series expansion of the exponential function in Eq. (|154p generates 
the perturbation expansion for the two-point Green function. In lowest order 
one obtains the non-interacting Green function Go whose inverse Gq 1 is 

G 1 (l,2) = (i 1 ^-m)5(l,2). 

The relation between a Green function and its inverse is 

/ d3G^ 1 (l,3)G(3,2) = / d3G(l,3)G^ 1 (3,2) = 5(1,2). 
Jc Jc 

Dyson-Schwinger equation is 

G-\l,2) = G 1 (l,2)-S(l,2), (156) 
which can be further written as 

G~ X G = 1 + SG 

GG l = 1 + GS, (157) 

whose explicit form is 

(«7 M ^-m)G(l,2) = 5(1, 2) + / d3S(l, 3)G(3, 2) (158) 

Jc 

G(l,2)(-i 7(J ^-m) = 5(1,2) + / d3G(l,3)£(3,2), (159) 

Jc 

After simplification by translating the contour integrals into normal ones we 
get Dyson-Schwinger equations 

rt 2 



(i 7A $* - m)G< (1, 2) = - / 2 d3S< (1, 3) [G> (3, 2) - G< (3, 2)] 

Jt 

- [ 1 d3 [E<(1,3) - S>(1,3)] G<(3,2)(160) 

Jtn 
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G < (l,2)(-i ll jj»-m) = - / 2 (i3G < (l,3) [E > (3, 2) -S<(3,2)] 

■ii.. 

- I 1 <& [G<(1,3)-G>(1,3)] S<(3,%61) 

where we have used the correlation parts of the self-energy and S < in 
analogy to Eq. (|155p . Making gradient expansion and Fourier transforma- 
tion, taking the difference of the above two Dyson-Schwinger equations, and 
finally taking traces for the resulting equations, we obtain Kadanoff-Baym 
equation 

d^Ti[^iG < (X,P)} = -Tr [iG > (-iS < ) - {-iT l > )iG < ] (X,P). (162) 

where X = (X\ + X2)/2 and P is the momentum variable conjugate to 
X. Note that our convention for the sign of selfenergy is implied in G^ 1 = 
Gf(°)-l - S, which means — iS in the Feynman rules. 

13.2 Kinetic equation 

We can derive the kinetic equation from Eq. Q162J) . The fermionic Green 
functions (here we mean the electron and the neutrino or their anti-particles) 
can be expressed in terms of distribution functions as follows 

iG<(X,P) = -^P^{f e -/At,P)S(PO+Ve-/u-p) 

"I 1 - fe+/u(t, -P)]<*(P0 +V>e-/v +P)} ■ ( 163 ) 

where we have shifted po by the chemical potential po — > po + fJ> e -/ u - In the 
same way, we can obtain G > as follows 

iG>(X,P) = ^P^{[l-f e - /v (t,PMPo + fi e -/v-p) 

-fe+/v(t,-p) S (jP0 + V>e-/v +P)} ■ (164) 

We assume that the distribution function / only depends on Xq = t, so that 
G(X,P) can be written as G(t,P). This approximation is well justified for 
neutrinos in the compact star with mean free paths comparable or larger 
than size of the star. 

We are interested in the neutrino cooling stage in the compact star. When 
neutrinos are produced they escape the star immediately. The neutrino 
density can not be built up in the star. So we can put the chemical potentials 
of the neutrino and anti-neutrino zero. Substituting Eq. 1|163|) and f)164p into 
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Figure 13: Self-energy in close-time-path formalism for the neutrino Urea 
process. 




Eq. (I165p . and carrying out the integral J °° dpo, we obtain the equation for 
the time variation rate of the neutrino distribution: 

B 1 f°° 

^/"(*>P) = 4^ I d PQ Tr[iG>(-iZ<) 

-(-iZ>)iG<](t,P), (165) 

which picks up the contribution corresponding to the term with support at 
Po = p in the left hand side. In the same way, we can extract the kinetic 
equation for the anti-neutrinos by carrying out the integral J° dpo 

§- t M^~P) = dp Tr[iG^(-iE<) 

-(-i£>)tG<](t,P), (166) 

which picks up the contribution corresponding to the term with support at 
Po = -P- 

In order to evaluate S > ' < , we need the Feynman diagram as shown in 
Fig. (fT3l) . The convention about the order of fermionic fields in the Feynman 
diagram is along the direction of momentum flow. From the diagram Fig. 
(fT3l) (a). the neutrino self-energy reads: 

-ffi<(p„) = -J ^- 4 (-irt)iGf(p e )ir» 

x[-m> u ,(P e - P u )]iDj^'(P e - P v )iD^'(P e - P v ) 

= 2j ^^(1 - 75)7^,^(1 - 75) 
x[-OI> (P e -P u )]-f e (t,p e )5(p' e0 -p e ), 

Pe 
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where the overall sign is due to the quark loop. Note that P e = (p' e0 ,p e ) with 
PeO = PeO+Me, then integral element is d 4 P e = dp' e0 d p e . We have defined the 
order of '>' and '<' in Il > ' < along the momentum direction. For example, 
in the first Feynman diagram Fig. (jT3]) (a). the W-boson polarization tensor 
is defined by from // to v' following the momentum flow, i.e. // as the first 
entry and v' as the second one. The first entry on the negative branch has a 
time later than the first entry on the close-time-path contour C. Therefore 
we have n^,. In the same way, Fig. (fT3l) (b) reads 



-«>(P„) = - J ^iY»iG>{P e ){-iYl) 



x[-m< u ,(p e - p u )\iD F /{p e - p v )w*S(p e - p v 



g\ r d 4 p P 



(27T) 



7 M (l-75)7 <T ^e, CT 7"(l-75) 



X[-0I< (P e - P„)]-[l " fe(t,Pe)W e0 - Pe 
Pe 

Inserting —iT^{P v ) and — zS<(P,) into Eq. (|165p . we get 
d . , s 1 Gl f dh 



dt 



x[l-Ut,p u )}f e (t,p e )[-iU> u (P e -P u )], (167) 



where P u = (p u ,Pv) and P e = (p e ,p e ). Note that the definition for P e has 
been changed from the original one. The leptonic tensor L^ v is defined by 

IT(P V , P e ) = Tr [y X P^(l ~ 75)7 <T P e ,.7 iy (l - Is) 

For the anti-neutrino, we still use P to denote its momentum, only the 
second term ~ (— iH~)iG~ in Eq. (I166|) contributes. Finally we obtain 

ot 4vr 2 J (2vr) 4 px7p e 

x[l - M-PuW ~ /e(Pe)]Hn< (P e - P v % (168) 

where Py = (—pu, Pv)- We can rewrite Eq. (I168|) by flipping the momentum 
sign: p F -> -Pf, 

|/,(*,p*) = -rffw^-- L "^- p '> 

ot Air 2 J (2vr) 4 p I7 p e 

x[l - / F ( PF )][1 - / e ( Pe )][-m< (P e - P F )], (169) 
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where Pp is changed to P v = {—pu, — Pf)- We will see that the negative sign 
in front of Eq. (|169l) will cancel the negative sign due to Pjt = (—pv, — Pf) 
in the matrix element. Thus the result for the anti-neutrino is the identical 
to the neutrino. 

Using the following relations between II> and n^, 

U<(K) = -2in B (k )lmIl R (K) 
U>(K) = -2i[l+n B (k )]ImU R (K), 

where n B (ko) = l/(e^ fe ° — 1) is the Bose-Einstein distribution, Eq. (|167p can 
be expressed in terms of Im.H R (P e — P u ) as 

^-Jv(t,Pv) = ]-G 2 F f ^--L^(P v ,P e )n B (p l/ -p e + (J , e ) 

Ot 47T J (2TT) /l p l/ p e 

X [1 - /„(*, Pu )]f e (t, Pe )ImU R u {P e - P v ), 

where we have used n B {ko)+n B {— ko) = — 1. Taking the momentum integral 
£^3 and noting that n v = f fv{t, P^), we obtain the number emissivity 
per volume for neutrinos: 

di Mt) = 2 Gf J (2.)3(2 7 t)3 2^2^ L {Pv,P*)nB(Pu-Pe + l±) 

X [1 - f v (t, P„)]fe(t, p e )ImnJ,(P e - P v ). (170) 

The emissivity for anti-neutrinos is the same as that of neutrinos. The energy 
emissivity can be obtained by including p v into the integrand in Eq. (I17Q|) . 

13.3 u-d-W vertex in Nambu-Gorkov basis 

Because we are interested in Urea processes, we have to write down the 
u-d-W interaction vertex in NG basis. The u-d-W vertex differs from the 
quark- gluon one in Nambu-Gorkov forms due to the flavor changing nature 
of the W-boson interaction. So we explicitly keep the flavor indices / for the 
quark fields: 

Mf) = C : f(f), ^ c {f) = ^ T {f)C 

where we assume / = 1,2 for u,d respectively. We need for the u-d-W 
coupling an additional flavor matrix t± to guarantee that couples to u 



68 



Figure 14: u-d-W vertices. The momenta of W-bosons flow into the vertices. 
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and d instead of the same flavor, where r± are defined by 
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Here we assume the spinor order of the vertex is d— W~ — u and u — W + — d, 
where the momenta of W-bosons flow into vertices, see Fig. (PT4|) . Then we 
have 



re 



2^2 sin 6 W 



y(i -75)t±, 



where 9\y is the Weinberg's angle. Here we have taken CKM matrix element 
for u and d quark V u d ~ 1. The Fermi weak coupling constant Gf is related 
to W by g§ = SM 4^e w - 

Then we can write down the u-d-W interaction term 



V>(/i)r^(/ 2 ) 



= ^(fi)CT^l(f 2 ) 

= ~ ? c (/ 2 )c T rf(7^ c (/ 1 ; 

= ^ c (/2)F±^c(/i), 

where we have defined = — CT^ T C. Therefore the u-d-W vertex in 
Nambu-Gorkov basis is given by 

ie 



± 



11 ±;NG 



2^2 sin 9 W 

ie 



2y/2 wiOw 
where we have defined 

7^(1-75) 



7 M (1 " 75)r± 

-7^(1 + 75K 

^NG T ±;NG, 



(171) 
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13.4 Polarization tensor for W-boson from quark loop in su- 
perconducting phase 

Now we can calculate the polarization tensor H^ u (Q) for the W-boson from 
the quark loop in imaginary time formalism for superconducting phases. 
Here Q = P e — P u . Then we make analytical extension to get retarded 
polarization tensor. The quark propagator with condensate can be written 
as Gij = GfjTn + GfjT22 where i,j are Nambu-Gorkov indices and t\\ = 
diag(l,0) and r 2 2 = diag(0, 1) in the flavor space of u and d quark. Gfj 
and Gfj are given by Eq. (|98l) . Let's look at one case where the left vertex 
is r_ ; 7VG while the right must be t + -ng, as illustrated corresponding to the 
fermion loop direction in Fig. (fT3l) . Then the polarization tensor of the 
W-boson in imaginary time formalism reads 

n""(Q) = JrW |^Tr WG , Ci/ , s [r^ G r„ ;iVG ^(p M )r^ G T +i v G g d (p d )] 

1 PuO J {l7T > 

= \ T E [ ^ T ^[r(i-i 5 )r-G^Pu) 

Puo J y ' 
xY(l-K)T + G d n (P d )] 

+ \ T H J ^|Tr cJiS [ 7 ^l + 75 )T + ^ 2 (P 11 ) 

PuO 

x 7 ^(l + 75 )r_^ 2 (P d )] 

PuO 

xya+^r^Pd)] 

-\ T Y, j ^Tr cJ , s [r(l + 75 )r + ^ 1 (P u ) 

PuO 

x 7 ^(l- 75 )r + ^ 2 (P d )]. (172) 

where P^ = P u + Q, and p u o = —i(2l + 1)ttT with integer I are Matsubara 
frequencies. The order in the trace is against the direction of fermion line in 
Fig. (jX3j) . The last two terms in Eq. (11721) proportional to the off-diagonal 
propagators are zero due to the presence of two identical r_ or r + leading to 
vanishing flavor traces. This a spectacular feature for single flavor pairing. 
Taking the flavor trace and Matsubara sum, and noting that the first and 
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second terms are identical, we obtain 



n^(Q) 



PuO 



d 3 p u 
(2vr) 3 



Tr C)S [ 7 ^(l - 75)^1(^)7^(1 - 75)^ (P d )] 



^Tr c>s [^(l - 75)70^" A P :"7^(1 - 75)7o^ rd Ap d ed ] 



n u (l-n d ) (l-n n )n rf 



+ 



-go + e« + £d 
(1 - n n )(l - n d ) 
-qo - e u + e d 



-qo — £u — £d 
n u n d 



-qo + e u -e d 



[1 - n F (e u ) - n F (e d )] 
n F (e u ) - n F (e d ){]>73) 



where V Tu is the projector corresponding to the excitation branch with index 
r u . A sum over the excitation branch r u is implied (the function inside 
the curly brackets depends on the excitation branch r u ). A p is the energy 
projector defined by A p = (1 — e7op • 7 )/2. n and 1 — n are Bogoliubov 
coefficients for quasi-hole and quasi-particle respectively. They are given by 



B^=n 
= 1 — n 



1 



2 + 2e 



(174) 



where £ = k — e/i and e = a/£ 2 + 2 - Note that the gap parameter eft depends 
on the excitation branch. In Eq. (|173|) we only keep particle excitations, i.e. 
with e u = e d = +. For the derivation of the Matsubara sum, see, for example, 
[871 [T4B]. 

Then we are in the position of obtaining the imaginary part of the re- 
tarded polarization tensor Imll by making analytic extension qo — ► qo — i5 
in Eq. (fT72l : 



ImlC(Q) = 7T 



d 3 p u 



TV C , S [7^(1 - 75)7o^ r "A- i( 7-(l - 75)70?™ A~J 



(2vr) 3 C ' SL ' v p " 
xn B 1 (-1o)Ba u B ad 5(-qo - a u e u + a d e d ) 
xn F (a u e u )[l - n F (a d e d )], 



(175) 



where a denotes the nature of the excitations: a = + for quasi-particles and 
a = — for quasi-holes. The factor n^{— qo) will cancel ns(p — Pi + /U e ) in 
Eq. (1170(1 . A sum over a u and a d is implied. 
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13.5 Evaluation of neutrino emissivity 

Substituting Eq. (|175p into Eq. (|170p . we get the number emissivity for 
neutrinos: 

d f d 3 p u d 3 p e d 3 p u 1 n 
—n v = Z7r / - — r; — — r= \M\ 

xB au B ad 5[p u - (p e - /i e ) - aie" + a 2 eif] 

xn F (p e - fi e )n F (a u e u )[l - n F (a d e d )]. (176) 

We have defined the matrix element 

\M\ 2 = G F p uPd Tt c ^(l- l5 h V r -A Pu 
X7"(l-75)7o^ rd A p jV(P„P e ) 
« G^Tr c , s [ 7 ^(l - 7 5)7o^ r "^,<x7 a 

X7"(l " 75)70^^71 p e) (177) 

where we have used the mass-shell condition for the u and d quark. The 
energy emissivity for the e-capture process is then 

9E U n _ f d 3 p u d 3 p e d 3 p u 1 2 



dt * J (2ir) 3 2p u (2n) 3 2p e (2tt) 3 2 Pu 2p~ d l 1 Pv 

xB au B ad 5\p u - (p e - fi e ) - a u e u + a d e d ] 
xn F (p e - n e )n F (a u e u )[l - n F (a d e d )]. (178) 

One can also obtain the time rate for the momentum emitted along any 
direction by replacing p v with p u cos 9, where 9 is the angle between the 
neutrino momentum and this direction. 

Before we evaluate Eq. f| 1 T8 [) . we analyze the phase space given by the 
energy conservation for quasi-particles, i.e. restricted by the delta function. 
If we don't take any interactions into account, for excitations near the Fermi 
surface, the energy conservation condition reads 

= p u ~ (p e - Me) - a«e« + a d e d 
~ Pv - a u (j) u + a d (j) d 
~ Pv, 

where we have used pi « [i% for % = u,d,e and p v 3> 4> u ^ d . We see in the 
above equation that in order to make the energy conserved, the neutrino 
energy must vanish, which means the phase space for the neutrino emission 
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Figure 15: Phase space for Urea process. 




is very much suppressed. So including the interactions for quasi-particles 
is necessary. The Fermi liquid correction to the non-interacting degenerate 
Fermi gas is a better description of interacting degenerate Fermi gas, which 
opens the phase space for neutrino emission. The Fermi liquid correction to 
the quark Fermi momentum is due to the quark-quark scattering via gluon 
exchange. The effect is that it reduces the Fermi momentum by an amount 
linearly proportional to the strong coupling constant as below: 

C F a s . 

PFi = (1 )VFi, 

where i = u,d. Actually the Fermi liquid corrections make it possible to 
form the triangle by pfu, PFd and pFe, as shown in Fig. (fl~5l) . We can then 
determine the angles B ud , B ue and 8 de by the triangle relation. For example, 
cos 6 u d is given by 



cos B ud 



Pld+Plu-Ple 
IPFuPFd 

(/xg + ^)(l-2«)-^ 

2/i d fi u (l - 2k) 
(2/4 + 2 Me /v)(l - 2k) + £{l - 2k) - 



2fj, d fj, u (l - 2k) 



H 2 e K 



Hd^u 

where we defined k = c F as . One can see that if one sets k = by turning 
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Figure 16: One way of simplifying the collision integral 



d 3 p u d 3 p d d 3 p u d 3 p 



dp v dp d dp u dQ d dB u 



momentum conservation 



b v dd t 



* — I — t 



u d0 u dp v dp d dp u 



freedom 2 



energy conservation 
freedom 1 



t 



collinear limit: d ~ 6 U 



numerical 



dp u dp d dp u dO d 



off the interaction, we have 6 ud = 0, which means the momenta of u and d 
are collinear. 

Another difficulty we have to tackle is the collision integral. Before inte- 
grating out energy-momentum conservation and making any simplifications, 
the integral is of 12 dimension. Carrying out the energy-momentum con- 
servation, we are still left with an 8-dimensional integral. For the normal 
quark matter, the rest integral can be carried out analytically. But for the 
color superconducting quark matter, one can only use numerical method to 
evaluate it. One way of simplifying the collision integral is illustrated in Fig. 
(US]). 

The neutrino emission and cooling rates for some of spin-1 phases are 
calculated in [1331 11341 1147] . The observational implication for the CSL 
phase has been studied in Ref. [14~8j 1149] . The basic ingredients in the 
calculation are the quasiparticle dispersion relations, containing the spin- 
one gap functions. We have studied in detail the effect of an isotropic gap 
function (CSL phase) as well as of anisotropic gap functions (planar, polar, 
A phases). The emissivity and the specific heat can be expressed in terms 
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Figure 17: Left panel: The functions G(tp,ip) of the neutrino emission con- 
tributions due to gapped modes in the CSL, planar, polar and A phases. 
Right panel: The functions K(ip) for four spin-one color superconductors. 




of the scaled gap parameter, cp = cp/T , 

G(ip u ,tpd) , 

(179) 

where e u and cy are the emissivity and heat capacity respectively, the nu- 
merical results for the functions G and K are presented in Fig. [17l 

In all four phases, also analytical expressions have been derived in the 
large cp limit (i.e., in the limit of small temperatures). In the case of an 
isotropic gap function (CSL phase), the well-known exponential suppression 
of the emissivity and the specific heat is observed. We find that anisotropic 
gaps give rise to different asymptotes in general. For example, the phases in 
which the gap has point nodes (i.e., polar and A phases) show a power-law 
instead of an exponential suppression at cp — ► oo. The actual form of the 
power-law depends on the behavior of the gap function in the vicinity of 
the nodes. While a linear behavior gives rise to a suppression ~ l/<p 2 , a 
quadratic behavior leads to ~ 1/cp. 

13.6 Anisotropic neutrino emissions in spin-1 transverse A 
phase 

As we mentioned that we can get the time rate of the momentum emitted 
along any direction due to the neutrino emission by replacing p u with p u cos 9 
in Eq. (|178j) . We can use the resulting formula to calculate the momentum 
changing rate for any color superconducting phases. Of particluar interest 
is the transverse A phase of spin-1 pairing. This phase shows an asymmetry 



457 „ 2 6 

——a s G F l [leV-u^d 



1 : 

3 + : 



Cy 



f=u,d 
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Figure 18: Effective gap of transverse A phase in neutrino processes, shaded 
area shows the magnitude of the gap. 




in neutrino emission along one spatial direction, which can be regarded as 
an example of parity violation in macroscopic stellar objects |150j . What 
makes the transverse A phase so special is that it shows a very spectacular 
helicity property in coupling with the neutrino emission. 
Let's look at the order parameter of the A phase: 



M = h 



7 I(k)+ 7 i(k) 



where 7_l = 7 — k7 • k. In the A phase, one direction is chosen to break 
the color and spin symmetry. This direction is assumed to be along the z- 
axis. 'V and '2' are two transverse directions to the z-axis. The excitation 
branches are characterized by eigenvalues A r of Ai^Ai. Here are excitation 
energies 



where A r are given by 



Ai, 2 = (l + |cos# fc | 2 ) 
A 3 = 0, 



where 6^ is the angle between the momentum and the z-axis. So we have 
spectral representation for Ai^Ai 

M^M = XiVi + \ 2 V2, 
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where the projectors are 



1 ± sgn(cos 9k)loTol ■ k 
= ijf [lTsgn(cos^)]i/ + (k) 
+^j|[l±sgn(cose fc )]ff-(k). 
Here the helicity projectors are given by: 

F±(k) = i [1± 75707 -k] 

where S = diag(<r, <r). We see that for the first excitation branch Ai, the 
mometum modes with the positive z-component must be negative in helicity, 
while the modes with negative z-component is positive in helicity. The sec- 
ond excitation branch just shows an opposite helicity correlation. Since the 
Urea processes involve exchanging W-bosons, only the particles with neg- 
ative helicity can participate. This fact can be seen by inserting H h with 
h = ± into Eq. (j!77|) where the Dirac trace becomes 

Quark Trace = Tr s [ 7 "(l - 75)70^ {Pu)A pu 

X7"(l-75)7o^ d (Pd)A p J- (180) 

We use 

S-pAp = -75A P 
F A "(p)A- = i(l-^ 75 )A- 

to rewrite the trace (I180p as 

Quark Trace ~ Tr s [7 M (l - 75)70(1 - /i u 7 5 )(p n )Ap t[ 
xt"(1- 7 5)7o(l -h d j 5 )A~ d ] 
= 0, for hu t d = +■ 

We see that components with positive helicity h U) d = + vanish. This selects 
upper part of the first branch and the lower part of the second branch, which 
gives an effective mode showing asymmetry along the z-direction, see Fig. 
(fT8l) . We know that the larger the gap is, the more suppression the neutrino 
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emission gets. So from the right graph of Fig. (fT8j) we see that the effective 
gap is larger in the upper half plane (positive z-component) than in the 
lower one (negative z-component). This means that the neutrino emission 
along the positive z-axis is more suppressed than in the negative z-axis. This 
asymmetry can be regarded as an example for parity violation in macroscopic 
scales. 
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